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A PFAFFIAN REPRESENTATION FOR FLAT ASEP 


JANOSCH ORTMANN, JEREMY QUASTEL, AND DANIEL REMENIK 


Abstract. We obtain a Fredholm Pfaffian formula for an appropriate generating func¬ 
tion of the height function of the asymmetric simple exclusion process starting from flat 
(periodic) initial data. Formal asymptotics lead to the GOE Tracy-Widom distribution. 
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1. Introduction 

In this article, we consider the one-dimensional asymmetric simple exclusion process 
(ASEP) with flat initial data, meaning that initially even sites are occupied and odd sites 
are empty. The particles then perform nearest neighbour asymmetric random walks, in 
continuous time, the only interaction being that jumps to already occupied sites are sup¬ 
pressed. ASEP is perhaps the most popular of the small class of special partially solvable 
discretizations of the Kardar-Parisi-Zhang (KPZ) equation. In recent years, a series of 
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breakthroughs on these models has led to exact formulas for one point distributions of the 
solutions of KPZ for very special initial data. The first formulas were for the narrow wedge 
and half-Brownian initial conditions [ACQ11; CQ13], which were obtained by taking limits 
of earlier formulas for ASEP with step and step-Bernoulli initial conditions due to Tracy 
and Widom [TW09a; TW09b]. The remaining key initial data, flat, half-flat and Brownian, 
are more challenging and had to wait several years ([BCFV14; OQR14], and the present 
work). Meanwhile there was a parallel effort in the physics community working directly 
from KPZ, using the method of replicas to derive rigorous formulas for the moments, and 
then writing divergent series for appropriate generating functions for the one point distri¬ 
butions of KPZ, which could then be formally manipulated into convergent Fredholm series 
[CDR10; DotlO; IS11; IS13], 

In this context, P. Le Doussal and P. Calabrese [LDC12] derived what in mathematics 
would be called conjectural expressions for the flat case, as Fredholm Pfaffians. To a 
certain extent, this work is an attempt to make sense of their remarkable, but highly 
non-rigorous computations. Our main results provide analogous formulas for flat ASEP. 
Formal asymptotics give the correct fluctuations in the large time limit, which in this 
case correspond to the GOE Tracy-Widom distribution [TW96], in a new form. In the 
weakly asymmetric limit one obtains versions of the formulas in [LDC12], However the 
technicalities of this, as well as those of providing a rigorous justification of the large time 
asymptotics, are so involved that we have left them for a future article. 

It should be noted that an earlier formula existed for the one point distribution of ASEP 
with flat initial data, due to E. Lee [LeelO]. However, this formula is similar to the first 
formulas of Tracy and Widom for the step case Tracy and Widom [TW08b] and are not in 
a suitable form for asymptotics, even at a formal level. In the step case, it already required 
an extraordinary argument [TW08a] to rewrite this as what is now referred to as a small 
contour formula, which is indeed suitable for asymptotics [TW09a]. In the flat case this 
method has not been successful so far, and thus we follow an alternative approach. 

In a prequel article [OQR14], we derived a formula for a generating function in the 
half-flat case, meaning that initially only even positive sites are occupied. Because such 
data is left-finite, meaning that there is a leftmost particle, the duality methods of [BCS14] 
apply, and with suitable guessing, inspired, in part, by the formulas of [LeelO], one is 
able to obtain an exact solution to the duality equations, which provide formulas for the 
exponential moments of the ASEP height function and can then be turned into a certain 
generating function. 

In this article we take a limit of the half-flat moment formula probing into the positive 
region, which reproduces the flat initial data. We follow the broad lines of [LDC12]. The 
asymptotics is extremely interesting, nonstandard, and unusually involved. One is forced 
to make a large deformation of contours, and the flat formula is an enormous sum of the 
residues passed as one makes this deformation. In this sum, a pairing structure arises 
which leads to the emergence of Pfaffians. 

Naturally, the remarkable formulas beg the question whether there is a more direct ap¬ 
proach, and some recent physics work have indicated possibilities in this direction [DNWBC14; 
CD14], 

Outline. The rest of the article is organized as follows. Our main results are discussed 
in Section 2. Section 3 contains the derivation of the flat ASEP moment formulas as 
a limit of our earlier half-flat formulas. Section 4 is devoted to unravelling the Pfaffian 
structure behind the moment formulas and rewriting them in a simpler way. The resulting 
formulas are then summed in Section 5 to form a generating function and showing how, 
in a certain case, they lead to a Fredholm Pfaffian. Section 6 contains a discussion of 
analogous results for the moments of the solution of the KPZ/SHE equation (or, more 
precisely, the delta Bose gas) with flat initial data. Appendix A presents a formal critical 
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point analysis supporting the conjecture that the ASEP height function fluctuations are 
described in the long time limit by the GOE Tracy-Widom distribution when the initial 
data is flat. As explained in Remark 2.2.vii, the rigorous asymptotics cannot be obtained 
directly by these methods and remains the subject for future work. Finally, Appendix B 
contains an overview of Pfaffians and Fredholm Pfaffians. 


2 . Main results 


2.1. Preliminaries. The one-dimensional asymmetric simple exclusion process (ASEP) 
with jump rates p € [0,1/2) to the right and q = 1 — p > p to the left is the continuous 
time Markov process with state space {0, 1} Z and generator 

= J2pv(x)( 1 - r](x + 1))(/(/? - U + l*+i) - f(v)) 

xgZ 

+ qr]{x + 1)(1 - r](x))(f(r] + l x - l^+i) - f(rj)), 


where l x (y) = 1 if y = x and 0 otherwise. Writing f){x) = 2g(x) — 1, we define the ASEP 
height function by 


h(t, x ) 


'ZN^it)+ Yf 0<y < x fi(t,y), x > 0, 
< 2A r g ux (f), x = 0, 

2N^(t) - Ex< y <o v(t,y), x < 


( 2 . 1 ) 


where Aq ux (t) is the net number of particles which crossed from site 1 to 0 up to time f, 
meaning that particle jumps 1 —>■ 0 are counted as +1 and jumps 0 —are counted as —1. 
For x ^ Z we define h{t, x) by linear interpolation. With these conventions h(t, x) is made 
of line segments with slopes ±1, and it evolves with its own Markovian dynamics, which 
flips V 4 A at rate q and A V at rate p. For q > p, the general trend is an upwards 
moving height function. 

An important parameter in the formulas will be 


T=p/qe[ 0,1], (2.2) 

We are primarily concerned with the flat , 

r^ at (x) = l xe2Z , 


and half-flat, 

Vo^ix) = lxg2Z >0 , (2.3) 

initial conditions (with the notation Z> n = {n+1, n+2,...} which will be used throughout). 
The superscripts flat and h-fl will be used for probabilities and expectations computed with 
respect to these initial conditions. 


In order to state our results we need to recall some definitions. Given a measure space 
(X,T,,p), the Fredholm Pfafflan (introduced in [RaiOO]) of an anti-symmetric 2 x 2-matrix 


kernel A'(Ai,A 2 ) = 


ATi(Ai,A2) A'i,2(Ai,A2) 
—A'i,2(A2,Ai) A~ 2,2 (Al 5A2) 


acting on L 2 (X ) is defined as the (formal) 


series 


Ff [j - K] mx) = £ [ MrfAi) • -MA*) Pf[/f(A„, Aj)]^,, 

fc >0 ■ Jx 

where, we recall, the Pfafflan of a 2k x 2k matrix A = {aij} is given by 


Pf[A] 


1 

Wk\ 


k 

Sgn(<r) a a ( 2 i~l),a( 2 i) ■ 
&&S2n i= 1 


Appendix B surveys basic facts about Fredholm Pfaffians. 
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We will also need various q- analogues of standard calculus functions. The q here is the 
q in quantum calculus [GR04]. Our parameter (2.2) came with the unfortuitous name r 
instead of q, since historically the left jump rate was called q. In an attempt to keep with 
standard practice, we will use q when discussing the definitions of the (/-deformed functions, 
but set q i—^ t within the computations. The standard q-numbers are given by 


[k]q 


Based on them one defines the q-factorial 



1 ~q' 


n q • — [1]<?[2]<jr • • • [n]q 


(g;g)n 

(1 - q y 


where the q-Pochhammer symbols are given by 


n —1 

(a; q)n = n ( X _ qka ) ■ 
k =0 

This last definition can be extended directly to n = oo when |</| < 1, yielding the infinite 
(/-Pochhammer symbols. The standard (/-exponentials are then defined as 

e q( z ) = ((i -q) z - q ) and Eq ^ = “ q ^ Z ' q ^°° ( 2 - 4 ) 

or (thanks to the g-Binomial Theorem) through their series representations 

_ k _ n ^k(k— 1 ) 7 k 

e i( z ) = J2in and E i( z ) =J2 k , ~ • ( 2 - 5 ) 

k> 0 q ' k> 0 q ' 


The first series converges absolutely only for \z\ < 1, but it extends analytically to < 0 
through the formula in (2.4), where it approximates exp(^) well as q 1. And, of course, 
E q (z) = e q(~ z )- I n [OQR14] we obtained a formula for E ll_fl [e T ((r h ^ ,x ^ 2 )J. 

In the flat case one is naturally led to use a different kind of (/-exponential. The symmetric 
q-numbers [NG94; FT95] are given by 


~ 0 fc/2 _ a -k /2 

^<7 = g i/2 _ q - 1/2 

and the symmetric q-factorial by 


with Oq! 


1. The symmetric q-exponential function is then defined as 


exp g (2) 


OO 


£ 



OO 


£ 


q\k{k-\) 

k q \ 


( 2 . 6 ) 


The function f(z) = exp g (z) is the unique solution of 

S -^ = m, /(o) = i 

Oq<, 

where the symmetric (/-difference operator is defined by 8 q f(z) = f{q l ^ 2 z) — f{q~ l ^ 2 z) 
(see for instance [GR04]). The standard (/-exponentials do satisfy analogous (/-difference 
equations, but, unlike them, the symmetric (/-exponential satisfies 


exp ? (z) = exp 9 -i(z), 

which is a natural symmetry in the flat problem (see Remark 2.2.v). 


(2.7) 
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2.2. Fredholm Pfaffian formula. Now we can state one of our main results. We need to 
introduce some functions (in all that follows we fix a parameter fi > 0): for s € C, A € M>o, 
y € C with \y\ = 1 and ( € C \ M>o, define 

^-/3, 2 +|s-i i — T s / 2 y (-r-*/ 2 ?/; r) 


^(s,X,y,C) = (-C Y 

and let, for ui € R, 


y 


1 + T s / 2 y (— r s / 2 y ; r) ( 


1+t _s / 2 i/ l + p/ 2 !/ 


727 


A ,y;C) = 


0 = 2 ^ / 
2^1 JiM 


ds e SUJ if(s,\,yX), 


( 2 . 8 ) 


(2.9) 


which should be regarded as the inverse double-sided Laplace transform of if (in s), see 
(5.7); additionally, if 2 , 21,22 € R and a, a\,a 2 € {—1,1}, let 

Fl _ V f-z) m (l - rl^^+2^-" 1 ^ 2 (^/VL (r 1+7n /?/ 2 ; r)^ 
11 ,yj J 1 J v 2 -1 (^x/ 2 ; ' r ) 00 { T /y 2 'i r )oo 

00 

^(* 1 , *25*1, * 2 ) = £ (-^) TO1 (-^) m2 (-aia 2 ) miAm2+1 sgn(a 2 m 2 -a 1 m 1 ) 

mi ,777-2 — 1 

r (|+/3)(m 2 +m 2 )-i(mi+m 2 ) 


(mi) T ! (m 2 ). 


w = -£ 


^2. r (|+/3)m 2 -im 


-Z 


m=l 


Theorem 2.1. For ( € C with |£| < r 1 / 4 , we /iaue 


E 


flat 


exp, 


.(^ t |M*>o ))1 =pf[j_ A' flat -C] 


( 2 . 10 ) 


( 2 . 11 ) 


1 L 2 (R>o) 

where the 2x2 matrix kernel Jf flat is given in terms of the functions defined in (2.8) and 
(2.10), with any (3 > 0, through (here Co,i is the unit circle centered at the origin) 


i t,C (Ai,A 2 ) = 


1 

+ 2 


dwidu 2 — -f dy$(uji, \i,y)$(u 2 , A 2 A)Fi(e 1 +UJ 2 \y) 

?. L 


dwi dw 2 V>(wi, Ai, ai)if(u2, X 2 , cr 2 )F 2 {e W1 ,e 01,02) 


o'i,o' 2 e{— 1 , 1 } 


< 2 U (Ai ; A 2 ) 

^ 2 a 2 t,C (A, ; A 2 ) 

Remark 2.2. 

(i) The symbol^ in the y integral appearing in the definition of 7L® at denotes that this is 
a principal value integral. This is because the y contour goes through the singularities 
of F\(z,y) at y = ±1 (see the discussion starting after (4.20)). 

(ii) The parameter (3 > 0 which appears in (2.8) and (2.10) is necessary to make (2.9) a 
convergent integral, but otherwise has no meaning. 

(iii) The three functions in (2.10) should be regarded as certain g-deformations of classical 

special functions. For example, F 3 (z) — 1 is a g-deformation of the exponential; in 
fact, F 3 (z) = 1 — exp T (— z) if (3 = 0. Similarly, F\ is related to the Bessel function 
Jo: 2 

F i ( z , y) = -T“T ( J o( 2 v^; F y 2 ) - 1), 
y z — 1 




duF 3 (e U )ip(u>-, Xi,a) 


( 2 . 12 ) 


o-e{-i,i} 


1 


= - sgn(A 2 - Ai). 
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(iv) 


(v) 


(vi) 


(vii) 


where Jo{z‘, t, y) is a 2-parameter (r € (0,1), y € C with \y\ = 1) deformation of the 
Bessel function given by J 0 {z-,T,y ) = Em=o( 1 - r ) 2 ™ r (§+ 2 fl™ 2 -™■ 

As y ^ ±1, J 0 (z-T,y) becomes J^\z;t) = Em=o T^ +2 / 3 ) m2 ~ m , which in 

the case ft = 0 can be regarded as a symmetric version of the Jackson g-Bessel 


function 


^\z-T) = Em=0 ( wWch ’ aS T 


1 as well, becomes the 


standard Bessel function Jq(z)). F2 is a bit more difficult to recognize, so we will not 
attempt it here, but in view of (146) in [LDC12] it should be possible to write it in 
terms of certain q-deformations of the Bessel and hyperbolic sine functions. 

The statistical symmetry of flat ASEP under the transformation (p,q) (q,p) to¬ 

gether with h i —y —h is respected by the formula because the symmetric ^-exponential 
satisfies (2.7). 

There are similar, although less appealing, formulas for generating functions of 
T2 h{o,x) d e f[ nec l i n terms of other q-deformations of the exponential function, see 
Theorem 5.1. These include formulas for the case of E T , but not e T . This last case 
does not seem to be accessible through our methods, because the 
are expected to grow like e cfe2 for some c > 0, and hence E flat 


moments of r 2 h (°’ x ) 
r (Cr 2 /i ( t ’°l)j can not 


be computed by summing moments. It is worth emphasizing here that the use of the 
function exp r (z) is intrinsic to the problem, and not just a convenient choice. The 
formulas for other generating functions are not Fredholm Pfaffians. 

The exp T generating function on the left hand side of (2.11) does not in general 
determine the distribution of h(t, 0). On the other hand, the formula in [LeelO] is for 
the distribution function of h(t, 0), but it is not obvious how to obtain (2.11) from 
Lee’s formula. 

Focusing on the left hand side of (2.11), it is tempting to believe that exp r (z) behaves 
sufficiently like exp(^) so that the key identity 


lim E flat 

t —>00 


exp{ 


e a(ft(t,0)-|i-P/ 3 r)| 


iiflat 



h(t, 0) — \t 
tV3 



(2.13) 


for a < 0 still holds with exp replaced by exp r . Indeed, this is the case with e T 
(this fact has been used succesfully in recent years to derive the asymptotics of some 
related models, although not for the half-flat or flat initial conditions, see for instance 
[BC14; BCF12; FV13]). 

For x > —(1 —r)^ 1 it is the case that exp r (^) looks quite like exp(z), and in fact it 
converges to it uniformly on [—a, 00 ) for any a > 0, as r / 1. However, exp T (z) has a 
largest real zero at xq ~ —(1 —r)^ 1 , and as x decreases below xo, it begins to oscillate 
wildly, with zeros at Xk € [xoq~^ k+1 ^ 2 , Xoq~ k ^ 2 ] for each k € Z>i - their precise 
disposition is unknown [NG94] - and reaching size approximately el log l x ^ / lo g( 1_r ) 
in between. Unfortunately, this genuinely precludes estimations cutting off this bad 
region. 

(viii) A formal steepest descent analysis shows that, setting £ = _ r - t / 4 +d' /3r / 2 ) the right 
hand side of (2.11) leads in the long time limit, as expected, to the GOE Tracy- 
Widom distribution. The limit is obtained in the form 


TgOe(t) = Pf[J - K r \ L 2([ 0jOO )) 


with 


K r ( A l5 A 2 ) 


tj(<9ai - d\ 2 )K A i(Xi +r, A 2 + r) Ai(Ai + r) 

2 Ai(A 2 + r) \ sgn(A 2 - A 2 )J ’ 


(2.14) 


(2.15) 


a formula for the Tracy-Widom GOE distribution which is essentially equivalent to 
one implicit in [LDC12], and which is also very similar to (but not quite the same 
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as) a formula appearing in [Fer04], Here the Airy kernel K A \ is defined as 

roo 

K Ai (\ 1 ,\ 2 )= di Ai(A!+0Ai(A 2 + 0- (2.16) 

J o 

A proof of (2.14) together with the formal t —)■ oo asymptotic analysis leading to this 
Fredholm Pfafhan is given in Appendix A. 

(ix) Although we have not included it, a similar formal asymptotic analysis in the weakly 
asymmetric limit leads to the Le Doussal-Calabrese Pfafhan formula for flat KPZ 
[LDC12] (see also Section 6 for related formulas). 


2.3. Moment formula and emergence of the Pfaffian structure. The sum in (2.6) 
includes a brutal cutoff by q n / 4 , so that the left hand side of (2.11) is practically a finite 
sum of moments. Thus, in a sense, the Pfafhan formula is mostly indicating a nice algebraic 
structure for a sum of moments. Therefore we state separately our moment formula. It 
involves a kernel A' flat which is related to, but not quite the same as the kernel R flat ’C 
appearing in (2.11). Since the meaning of the terms in A' flat only becomes apparent in the 
computations, we do not repeat the detailed formula here. 


Theorem 2.3. Let i7 flat be the the 2x2 matrix kernel given by (4.24) and (4.26). Then 
for any m € Z>o we have 


E 


flat 


T \mh(t,0) 


= m T W 


171 / t\k 

.-Vrt" 1 ) 


k =0 


k\ 


oo ,. 

V / dAPf[A' flat (A Q ,A b ;m a ,m fe )]^ =1 . 

™.-i 4(K >0 )fc 


mi H-1 -mk=m 


(2.17) 


The route to these moment formulas for the hat initial data is by taking appropriate 
limits of formulas for the half-hat initial case (2.3). The limit we are interested in consists 
in starting ASEP with the shifted half-hat initial condition 770 ( 2 /) = l- y e 2 Z>_ x , considering 
the variable h(t,0) and computing the limit when x —>• 00 . More precisely, we will obtain 
moment formulas for hat ASEP through the identity 


jgflat 

T \mh(t, 0) 

= lim E 2Z> ~* 

T ±mh(t, 0) 



x—foo 

. 


(2.18) 


for m> 0, where the superscript on the right hand side simply refers to the initial condition 
specihed above. Introduce now the random variables 


X 

N x {t) = Y vt(y), 

y=-oo 


which are of course finite if and only if the initial condition is left-finite. It is not hard to 
check that when all particles start to the right of the origin one has N^ ux (t) = No(t), and 
hence from (2.1) we get 


h(t,x ) = 2 N x (t) — x (2.19) 

in the half-hat case. On the other hand a simple coupling argument shows that h(t, 0), 
with initial condition 770 ( 2 /) = l y ^ 2 Z > - x has the same distribution as the shifted observable 
h(t,2x) with initial condition 770 ( 2 /) = l ye 2Z >0 - Hence, using (2.18) and (2.19) we deduce 


that 


•jpiflat 


T \mh(t, 0) 


lim E h ’ fl 

x—too 


T m(N 2 x{t)-x) 


( 2 . 20 ) 


In Theorem 1.3 of [OQR14] we obtained the following formula for the moments appearing 
on the right hand side of (2.20): 


Theorem 2.4. Let m > 0. Then 


E 


h-fl 


r mN x (t) 




k =0 


= m- 


( 2 . 21 ) 
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kl 


L 


dw det 


niH-i-nfc=m 


rik 

O 0,7-^ 


-1 


W a T' la — ICft 


- a,b=l 


x n f(^a;n a )g p (u; a ;n a ) f)i(w a , w&; n a , n b ), ( 2 . 22 ) 

a=l l<a< 6 </c 


where C 0t ~ v is a circle of radius t 77 centered at the origin, with r) € (0,1/4), C k )r _ ri 
denotes the product of k copies of C 0 t -h , and where 


f (w;n) 


(1 _ T^^i-pAih 



x—1 


and 0 P and f)i are given in (2.23). 


For simplicity, throughout the rest of the paper we will omit the bound on the indices 
in products such as 1 < a < k and 1 < a < b < k when no confusion can arise and 
the factors involved in the products are defined in terms of a collection of k variables. A 
similar convention will sometimes be used for sums. Additionally, we will continue using 
the notation C k for the product of k copies of a given contour C in the complex plane. 

One of the most interesting parts of the story is the computation of the limit (2.20), 
which involves taking the limit x —>• oo of the right hand side of ( 2 . 22 ) with x replaced 
by 2x and after multiplying by T~ mx . We will start by rewriting this formula in a slightly 
different way which will be better adapted for the calculation of this limit. Introduce the 
following functions: 






h( w : n ) = 


_ / 1 +r n w ^—-kn 


l+«j 


T 


Q p (w;n) = 
flu (w a ;n a ) = 


(~w,r) o', ( 


T 2 n W 2 \T 


)« 


(-t*w,t) 00 (r n u; 2 ; t)^ ’ 

(-Wa,r) oo (r 2na u; 2 ; r) c 




(■ -T na w a ] t)^ {T l + n -w 2 -, r)^ ’ 
(w\W 2 ] r)oo (T ni+n2 w 1 w 2 ] t) c 


^(wajWb -,n a ,n b ) = 


( T ni w\w 2 ; t)^ ( t U 2 w\w 2 ; r ) 0 
[w a T na - W b T nh )(w b - Wa) 


(' W a T na - W b )(w b T nb - W a ) 


2x-l 


(2.23) 


g u will not appear in the coming formula (2.26), but will appear later in Theorem 2.5. The 
subscripts in g p and g u stand for “paired” and “unpaired”, a terminology which we will 
explain shortly and which will become clear in Section 3. Using these definitions in (2.22) 
(with x replaced by 2x, and multiplied by r~ mx ) and expanding the determinant in that 
formula using the Cauchy determinant formula 


we obtain 


det 


1 


x a - Vb 


a,b= 1 


Ua <b (Xg ~ X b ){y b - y a ) 
n a,b( X a - Vb) 


E- 


h-fl 


r m(N 2x (t)—x) 


= m T \ 


2X ) 


k =0 


(2.24) 


(2.25) 
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with 


^ ^ r _ r 

^k%(^ 2x ) = Y / 27 ri) A; / k d wYlh(w a -n a )f 2 (w a -,n a )Q p {w a -,n a ) 

m,...,n k >l ' ; C o,T-v a 

ni~\ -frifc=m 

X JJ t)i(w a ,w b ;n a ,n b )l) 2 (w a ,w b ;n a ,n b ). (2.26) 

a<b 


We have added a tilde in 2x) to indicate the fact that we have multiplied by r mx 

in (2.25) (in other words, v^f n {t. 2 x) = r~ mx n^(t,2x) in view of (2.21)). Note that in 
(2.26) we have written r~ mx as [7„ r~ naX . 

Note that the only factors in the integrand on the right hand side of (2.26) that depend 
on x are those of the form f 2 (w a ',n a ). An easy computation shows that the base in this 
power has modulus strictly less than 1 if and only if |rc a | > r _Tl “/ 2 . This suggests that 
for each a = 1,..., k we should deform the corresponding contour C p T - v to some contour 
lying just outside the ball of radius r” n “/ 2 . 

As we perform this deformation we will cross many poles. The residue calculus associated 
to this deformation is quite complicated, and is explained in detail in Sections 3.1 and 3.2. 
The result has two properties which turn out to be crucial for the sequel. First, there 
are two types of poles that are crossed (see the list in page 11): the ones coming from 
0 p (u; a ,n a ), which occur at w a = ±T~ n “/ 2 , and the ones coming from fq (w a ,Wb',n a ,nb), 
which occur at w b = T~ Ua /w a whenever n a = rib■ We will refer to these two types of poles 
respectively as unpaired poles and paired poles. This pairing structure is the key to the 
emergence of the Pfafhan. 

The second crucial property is that the residues of the unpaired and paired poles are such 
that x disappears from the corresponding factors (see (3.2) and (3.4)). As a consequence 
of this, after deforming the contours x only appears in factors of the form f 2 {w a \n a ) with 
w a living in the deformed contour, which are such that |f 2 (^a!^a)| —> 0 as x —» oo. This 
will allow us to derive an exact formula for the limit of 2 x), which after further 

simplification and a change of variables which turns the contours into circles of radius 1 
(see Section 3.3), leads to the result that follows. To state it we define the functions 


(1 — T ni WlW 2 ){l — T n2 WlW 2 ) 

e(wi,w 2 ,n 1 ,n 2 ) - ^ _ WlW2 )(\ _ T ni+n2 wiw 2 ) 
b(wi,w 2 ;ni,n 2 ) = t]i(w 1 ,w 2 ,n 1 ,n 2 )i) 2 {w 1 ,w 2 -n ll n 2 ). 


(2.27) 


Furthermore, given generic functions a(z,n) and b(zi, z 2 ; n 2 ), we define modified func¬ 

tions a and b through 

a(z,n) = a(r~ 2 n z,n) and b(zi, z 2 \ n l5 n 2 ) = b(r~ 2 ni zi, T~^ n2 z 2l n\, n 2 ). (2.28) 

In the formula that follows this modifier will be applied to all the functions defined in (2.23) 
and (2.27) (the modification arises through the change of variables performed at the end 
of Section 3.3, see (3.18) and (3.19)). 


Theorem 2.5. Let 


= lim u k%{t, 2x )- 

nr .—VfYi ’ 
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Then, with the notation introduced in (2.23), (2.27) and (2.28), and letting Co,i he a circle 
of radius 1 centred at the origin, 

1 


■r 


“ (*) = E 


E E E 


l 2 fc P k ! /—j "Ea n S+ 2 E a nl=m 

k u ,k p >0 P cri,...,crj. u G{—1,1} 1 n^,...,Ti^ >1 

fc u +2fcp=fc P 

k it 

[ kp d ™ P II(fi0u)(CTa;n a )]^(figp)(zP;nP)(fi0 p )(^;nP) 
l 7rl / "'Co, 1 ! a=l a=l 


k u 


X 


n^ ri 3 ” s n 


( 1)T V-*"S("S+i)(T-"S)ii 2( *P ,4;nP,nP) 


a=l a=l 

X h( 

l<a<6</c u 






n 

a<k u , b<kp 


t <n >>l(a a , z p ; <, n£)c(<r a , <, n£) 

b 


(2.29) 




In view of (2.20) and (2.25), this provides a first formula for °)/ 2 ], which will 

be later turned into (2.17). A minor issue with this formula is that the numerator and 
denominator in the factors of the form f}(cr a , cry nJJ, n'f) vanish when a a a b = 1 and — nj) 
is even, in which case the factor is to be interpreted as in (3.5) below (see the discussion 
preceding that definition). 

This formula can be regarded as a microscopic (and rigorous) version of formulas (107)- 
(108) of [LDC12], which can be recovered in the weakly asymmetric limit (see also Section 
6 ). 

To get an idea of how the pairing structure arising from the contour deformation leads 
to a Pfaffian formula for let us consider first the product 


cOM; ^, P ; nl,nl)i(z a p , z p ; n,P, n p )l(z a p , z p ; n p ,n p ). (2.30) 


l<a<6<fc p 


It turns out that this product has the structure of a Schur Pfaffian. To see this, we recall 
first Schur’s Pfaffian identity. 


Letting y a 


n 

l<a<b<2k 


X a ~ X b 


X a + X b 


Pf 


Xg ~ X b 
X a + X b 


2k 

a, 6=1 


Xg — 1 
Xa ~\~1 


we obtain the identity 


Pf 


Vb ~ Va 

MaVb ~ 1 


2k 

a, 6=1 


n 

l<a<6<2fc 


Vb ~ Va 
VaVb ~ 1’ 


(2.31) 


(2.32) 


Although the four factors inside the product in (2.30) are not of the form (y b —y a ) / (y a y b — l), 
this expression appears after rearranging the factors in the numerator and denominator. 
In fact, letting 

V 2 a -1 = z p and y 2a = z a p 

for a = 1,..., k p , one checks that the product in (2.30) is given exactly by n a <& yay b -i ’ 
and thus it equals the Pfaffian in (2.32). 

In Section 4 we will show how the last two lines of (2.29) (which include in particular 
the product in (2.30)) can be turned into the product of two Pfaffians, times the products 
of the factors of the form T ~ na - n <> and r~~ aU b . We will then show how the two types of 
variables (paired and unpaired) can be put on the same footing, and finally how the whole 
integrand in (2.29) can (almost) be written as single Pfaffian. The result is (2.17). From 
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this we can form a generating function. This is done in Theorem 5.2, a particular case of 
which becomes, after further manipulations, Theorem 2.1. 

The Pfaffian structure which we have described was first discovered in the context of the 
delta Bose gas with flat initial data in [LDC12], using a rough singularity analysis. They 
started from the divergent series for the generating function, while we start from moment 
formulas (4.18) which are more reminiscent of the Borodin-Corwin approach [BC14]. In 
hindsight, the two give rather similar results. However, the existence of nice exact expres¬ 
sions for the flat case at the ASEP level is a surprise, and it is far from clear that it is 
replicated in other solvable models such as q-TASEP, the log-Gamma polymer, and the 
O’Connell-Yor semi-discrete polymer. 

We next describe in detail the core of the argument, which is the residue summation. 


3. Moment formula: contour deformation 

And though the poles were rather small 
They had to count them all 

The goal of this section is to prove Theorem 2.5. We are trying to compute the limit as 
x —» oo of where is defined in (2.26). As we noted in the introduction, the only 
factor in the integrand on the right hand side of (2.26) that depends on x is f 2 (w a ,n a ), and 
the base in this power is strictly bounded by 1 in modulus whenever |u> a | > T ~ n o./ 2 _ Hence 
in order to compute the limit, we will deform the w a contours to some contour which has 
this property, which for simplicity we choose to be the contour Cq “ given by a circle of 
radius i-~ n °-/ 2 ~ T i centered at 0. 


3.1. Poles of the integrand. In order to write down the result of the contour deformation 
we need to determine first what poles are crossed as we deform each of the contours. It 
will be useful to list first the possible singularities of the integrand: 


(a) fi (w a ,n a ) has a simple pole at w a = 0, but this point is inside the original contour 
C 0 t - v . It also has essential singularities at w a = — 1 and w a = — T~ n “, but the first 
one is contained inside C 0r -ri while the second is outside Cq 0 . So this factor will not 
contribute any poles. 

(b) f 2 (w a ,n a ) has simple poles at w a = —1 or w a = —T~ na (depending on the sign of 
2x — 1), but as in (a) they are never in the deformation region. 

(c) g p (w a ,n a ) has simple poles at w a = —r~ 1la ~ t for any £ € Z>o, but all these points lie 
outside CgY It also has simple poles at w a = ± T “ n “/ 2 ^ f / 2 for any £ € Z>o- These 
points are inside the deformation region only when £ = 0, and we have, for a = ±1, 

Res i 9 p {w k ,n k ) =-±<JT-^ nk g u (aT-^ nk ,n k ). (3.1) 

Wj c =<7T—'2 ri k 


Observe also that, crucially, for a = ±1, 


f 2 (o-r 2 nk 



1 +CTT n k/ 2 ~bn k 
1 +<JT- n k/ 2 


2x—l 


= a. 


(3.2) 


(d) The factors on the denominator of hi (w a ,Wb',n a ,nb) vanish when w a Wb = T~ na ~ ia 
with i a € Z>o and w a Wb = with -4 € Z>o- Observe that whenever one of the 

two factors vanishes, the factor (yj a u’b',T) 00 in the numerator of hi (w a ,Wb',n a ,rib) is 
going to vanish as well. Thus there can only be a (simple) pole when both factors on 
the denominator vanish at the same time, i.e., when Wb = w ~ 1 r~ na ~ f - a = w~ 1 r~ nb ^ £b 
for some l a , 4 € Z>o- Whether these points lie inside the deformation region or not 
will depend on which stage of the deformation we are at: 

(d.i) Assume first that w a € C 0t - v and we are deforming the Wb contour. The possible 
singularities occur at Wb = tc“ 1 r _nf,_£6 for some 4 € Z>q. This point has modulus 
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T v~ n b-tb > T ~ 2 n b~v (recall r/ G (0,1/4)) and thus it is outside the deformation 
region. As a consequence, this factor does not contribute a pole in this case, 
(d.ii) Now assume that w a G Cq“ as we deform the w b contour. There is a singularity 
when Wb = w\ = w 2 with w\ = and w 2 = Observe 

that lu) 1 ! = T n o./2+v~ n b-h w hich is smaller than T~ nb ^ 2 ^ ri if and only if £\ < 
(n a — rib )/2 + 2 rj. Similarly we have |u) 2 | = T ~ n o./^+ri-h which i s smaller than 
T ~ n b/ 2 ~ r i if a nd only if £ 2 < (n b — n a )/2 + 2 r). Remembering that we also need 
to have £i,£ 2 > 0 these conditions on £±,£2 (and the fact that r/ G (0,1/4)) imply 
that n a = rib , and thus also that £1 = £2 = 0. Hence we have a simple pole at 
Wb = u;“ 1 T -na when n a = rib, whose residue is computed to be 


Res bi{w a ,w b -n a ,n a ) 

w b =w a 1 T~ na 


(—l) na+1 r — n a (n a +l)/2 ( -n a 
W a ' 


I)- 


(3.3) 


Note that, crucially, in the case n a = rib we have 


h{'w b 1 T 2 nb ,n b )f 2 (iUb,n b ) 


( T+W^T 2 T + T n bw b -n b 
l r+u;“ 1 r 2 -"i> r+w b 


(3.4) 


(d.iii) The third possibility of interest is when w a = or~ na ! 2 for a = ±1 (which happens 
when a residue in w a has already been evaluated as one of the poles from (c)). 
Exactly as in (d.ii), a singularity may only arise when n a = n b at w b = or~ na l 2 . 
But in this case the factor f) 2 (^6, ctt”"" 1 / 2 ; n a , n a ) has a double zero at this point, 
and hence the double zero in the denominator is canceled, which means that there 
is no pole in this case. 

(e) Finally, the factors on the denominator of t) 2 (w a ,w b ;na,nb) vanish when w b = w a T na 
and w b = w a T~ nb . The first point is inside C 0jT -r, whenever w a is in C 0 r ->7 or C^ a . 
The second point is outside C// b whenever w a is in C () T - V or Cq 11 . So this factor will 
not contribute any poles. 

Finally, consider the product (fjifoX^a, (T b T~ nb ^ 2 ] n a , n b ) at w a = o a T ~ na / 2 for a a ,a b = 
±1, which occurs when both w a and w b are evaluated as residues coming from the poles 
in (c), with w b having been evaluated first. When n a — n b is even and a a a b = 1, both the 
numerator and the denominator of the factor coming from [)i vanish. Since the evaluation 
at these points comes from a residue computation (in w a ), the product should be computed 
as 

(f)if)2 )((JaT~^ na ,abT~2 nb -n a ,n b ) := lim (f)ifl2 ){w,a b T~^ nb \n a ,n b ). (3.5) 

W — yO’a'T na ' 2 


The right hand side is computed in the next lemma. We stress that, in the formulas that 
will follow, we will use the extension of fqfo defined in (3.5) mostly without reference. 


Lemma 3.1. For n a ,n b G Z>o and cr a ,cr b = ±1 we have 


lim (fhf) 2 )(w,a b T 2 nb \n a ,n b ) 

W —(jT na 


{-a a a b ) n ^nb T -\nan b 


\a b A n b-a a A n °-\ 

1 -a a a b A^+ n b) ' 


Proof. We rewrite the limit as 

lim(fht)2 )((r a T~^ na w,a b T~^ nb ;n a ,n b ). 

w—ll 

A simple computation (see the proof of Lemma 3.5 for a similar one) leads to 

2 na ~ 2 nb y) 
T t-\n a +\n bijj 


(f)if) 2 )(v a r z na w,(j b T 2 nb \n a ,n b ) = \ 


1 - 


1=1 


1 - 


(3.6) 
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where w = a a abW. This product equals 


pj- 1 - T l ~ l+ h n --h nb ij, 


na w - r ^-i-5 n “+5 n & 


1 — 2 na +2 nb UJ 1 — 2 na+ 2 nb W 


n 


—T 


\-Z+\n a -\n h 


1 1 W — T 2 


b 9 


n"= 2 ( w - r 


•—1— hria+brib 


1 — T2 


inc+ln^ _ ^-ina+lnj^ 


The last ratio equals n(=i 1 crarTb f’ 1 T —r-——, and each factor in this product tends to 

\—T 1 2 na + '2 nb craCr b w 

cr a &b as w —> 1 unless (-—\n a +^nb = 0, in which case it tends to —1. Therefore the whole ra¬ 
tio goes tO (<7 a <T)) n “~ 2 (— l) lna ~ n b £27 ‘>0 (a a (Jl ) ) 1 ~ lrla ~ rlb£27 ‘ >0 = (<7 a (7fe) n “~ 1 ( —(T a lTh) lria_n () 62Z >0, 
and thus the limit in (3.6) equals 


_ (— a a ,) n a + 1 na-n b £2Z >0 -^ T ^na-^n a n b CT a C T b ~T? nb 

^ “ ' l-r^+^ b a a a b 


= 


— T \ n a I 


. l T Wb_ T W 

(_an«Arif, — j n a n b \ 

l — ri^a + ^b 

~-ina.n b T ^i>+ri n “ 

T Z -T—— T 

l_|_ r 2 na + ^ n b 


if a a (Jb = 1 
if a a a b = -1, 


which gives the claimed formula. 


□ 


3.2. Contour deformation and computation of the limit. As we mentioned, we will 
move the contours one by one, starting from Wk, then Wk—i, and so on. From the above 
discussion, there are two types of poles which are crossed. First, as we expand the w a 
contour we will cross singularities at w a = Tr - ”" 1 / 2 as in (c) above. We will refer to these 
as unpaired poles. The n variables associated to these poles will be called unpaired variables , 
and will be denoted with a superscript u (note that unpaired poles have no w variables 
associated to them, since a residue has been computed). 

The second type of pole appears at w a = w^ l r~ nb (for a < b) when Wb € Cq 6 and 
n a = rib- We refer to these as paired poles. To each paired pole corresponds a pair of w 
variables and a pair of n variables: Wb, which is integrated over Cq 6 , and w a , which is set 
to w a = Wfr 1 T~ nb , and similarly and n a , which are set to be equal. We call these paired 
variables, and denote them with a superscript p. 

Finally, in addition to the paired and unpaired variables, there are free variables , denoted 
with a superscript f, corresponding to the w variables (and associated n variables) which, 
after the deformation, end up integrated on a Cq contour (note that these have no paired 
variable associated to them). 

The result of deforming all the w a contours involves a sum over all possible ways in 
which the k original variables can be grouped into k u unpaired variables, k p pairs of paired 
variables, and fcf free variables, which of course means that k u + 2 k p + k{ = k. Given such 
a grouping, the original Wa s and n a ’s will be split into strings of variables: n u for the 
unpaired ones (of size k u ), n p and w p for the paired ones (of size k p ), and n f and for 
the free ones (of size kf). It is important to note that in the case of the paired variables, 
the variables forming these strings will appear twice in the formulas that will follow, once 
associated to w p and once associated to the “paired” string w~ p defined by 

w 7 = T ~ n "/ W a for a = 1,..., k p . (3.7) 

These variables should be regarded as “virtual”, in the sense that they are not being directly 
summed or integrated, but are instead defined directly in terms of other variables. On the 
other hand, in the case of the unpaired variables there is an additional string a of ±1 
variables, which will indicate whether the unpaired pole is evaluated at r _Tl “/ 2 or —r~ na ^ 2 . 
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It will be convenient to introduce another string of virtual variables w u , defined as follows: 

for a = 1,... , k u . (3-8) 


11 — —77 U 

= (T a T 2 n a 


Before writing down the result of the contour deformation we need to introduce some 
additional notation which will help us express some products arising from the pairing 
structure which emerges from the contour deformation. Assume we are given a collection 
of pairs of strings of unpaired variables (u; u , n u ), paired variables (uP,n p ) and free variables 
(«J f ,ra f ). Given a function g of a single w and a single n we define 

fcu kp fcf 

V u,p,t g(w u ,w p , ?u f ; n u , n p , n f ) = g(w^,na)]^ g{w p ,n p )g{w a p ,n p )]^ g{w[,nl). (3.9) 

a= 1 a= 1 a =1 

That is, we multiply over all possible variables, accounting for the extra set of paired ones. 
Similarly, given a function h of a pair W\,W 2 and a pair ni,ri 2 we define 

V u ’ p ’ i (w u ,w p ,w i ;n u ,n p ,ri i ) = JJ <, n£) JJ h(r^, n^, tt|) 

l<a<b<k u l<a<b<kf 

X n h ( w a,wl-,n p ,n p ) h ( wP ,w b p -, n p , n p )h{w a p ,w p ; n p , n p )h(w a p , w b p ; n p ,n p ) 

l<a<b<k p 

x II h (K, w b'’ n a, n bM W a’ W 7'’ n ai n b) II h ( w ai w h n a, n b) 

a<k u ,b<k p a<k u ,b<k{ 

x II Kwl,wl\n p ,n[)h(w~ p ,w[]n p ,n\). 

a<k p , b<kf 

(3.10) 

Here we multiply over all possible pairs of different variables, accounting again for the 
extra set of paired variables, but omitting factors with variables from the same pair (that 
is, omitting factors of the form h(w p ,w a p ;n p ,na )). 

We extend the notations (3.9) and (3.10) as follows: if any group of variables is omitted 
from the superscript (and the argument), it simply means that those variables are omitted 
from the products. Thus, for example, 

k 

V u g(w u -, n u ) = f[ O (3-11) 

a= 1 

while 

V n ’ p h(w u ,w p -,n n ,n p ) 

= II II 

l<a<b<k u a<kuib<kp 

X II h ( W %’ W b 5 n l)K W ai W 7’ n l)K w 7i W b ; n E> n b) h ( W a P , W 7' O' 

l<a<b<kp 

(3.12) 

We will slightly abuse this notation by omitting the superscript in V (which in any case is 
implicit in the superscripts of the argument). For convenience we will also write 

^(•;-) = ^i(-;-)^ 2 (v)- (3-13) 

Finally, we introduce the index set 

Ktote = {(<?,n\n 2 ,n 3 ) € {-l,l} fcl x Z* x Z^ x 
k\ k2 ks 

Y2 + 2 n l + Y2 n a = m, v a n l a + cr b nl V 1 < a < b < ki } (3.14) 


a=1 


a= 1 


a= 1 
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and the notation Cq = C p 1 x • • • x Cg fe for n of size k. 

We are ready now to state the result of the complete contour deformation. 

Proposition 3.2. For every k,m € Z>o with k < m we have 




f n (t,2x) = J2 


1 


3m fl (^u^p,fcf) 


A^U 5 Alp 5 Alf ^ 0 
Aiu "l - 2 Aip "(“/cf ^ = k 


k u \ 2 k p kJ.kfl 


P) n 'Wi 


where z/^ was defined in (2.26) and 


^(k u ,k p ,kf) = 


(o , ,n u ,n p ,ra f )eA™ fcp fcf 


( 2 vri) fc p +fc f J c * p xcf 


L 


dw p dw { Pfi(w u , w p , w f ; n u , n p , n f ) 


x Pf 2 (w f ; ra f )P 0 u(^ u ; n u ) 7 :, 0 p(w p , w f ; n p , n f )"Pfi(w u , w p , w f ; n u , n p , n f ) 




a=l a=l 


(3.15) 


The proposition follows from a slightly more general result which we prove next. Define 
the annulus Rfi = {u> € C: r _?7 < |w| < r -n / 2-r? } and the set D = IIn>i -^o x ( n }- We 
assume as given functions q i: D —>• C and i/j: D x D —>• C such that n ) is analytic in 

w € i?g and tf(wi,W 2 ',n\,n 2 ) is analytic in each variable wi € Rfifi and W 2 € i?g 2 (we are 
slightly abusing notation here to write 'ip(wi,W 2 ',ni,n 2 ) instead of n\), (w 2 ,ri 2 )))- 

Proposition 3.3. Fix k,m € Z>i with k < m and functions as <f> and if as above. Then 

e j^k L d ™ n </>(w a , ^a)flp(^a, n a ) JJ ^{w a ,w b \ n a , n&)f)i(w a , n a , n&) 

ni,...,nfc>l ^ ' ^o,r — ^ a 

niH-} -n^—nn 


a<b 


E 


A:! 


E 


1 


ri \kp-\-kf 


f dw p f 

Jcf p Jcz 


dw 


k u \2 k v kJ-kfl 4^ (27ri) 

fc u ,fc p ,fc f >0 u p (a,n“,nP,S f )eA^ kkc 
k u -\-2k p +kf=k u ’ p ’ f 

x Vfi(w u ,w p ,w i ] n u , n p , n i )Vfi>(w u , w p , w { ; n u , n p , n f ) 
x P 0 u (w u ; n u )P 0 p (w p , w f ; n p , n f )Pl)i(w u , w p , w f ; n u , n p , n f ) 


x II ^ Il(- 1 )'"“' r 

a =1 a=l 


nS T -ng(ng +1 )/2 (1 _ <P; „P, „P). 


Proof. The proof is by induction in k (for fixed m). First we check the case k = 1. Note 
that Vt)i(w’,n) = 1 in this case (because the product is empty), and thus the only poles 
crossed as we deform are the ones coming from 0 p (wi,ni) as in (c) in page 11. On the 
other hand, the condition k u + 2 £; p + kf = 1 implies that the sum on the right hand side has 
two terms: one with k u = 1 and the other two set to 0 , and one with k u = 1 and the other 
two set to 0. Since f>(wi,ni) is analytic on the region of interest, one checks directly that 
the first term corresponds to the computation of the residues at the two poles (keeping in 
mind that the poles are passed from the outside, so we get an additional minus sign), while 
the second term corresponds to the integral on the deformed contour. 

We assume now that the formula holds for a given k and rewrite the integrand in the 
case k + 1 , for fixed Wk+i € C 0 t -» 7 , as 

k 

(4>Qp)(wk + i,nk + i)Y\_(4>9p)( w a,na) (4b i){wa,w b ;n a ,n b )(w a ,w b ;n a ,n b ) 

a= 1 l<a<fe<Ai 
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with 

4>(w a ,n a ) = (j>(wa,n a )(^bi)(wa,Wk+i;na,n k+ i). 

By (d.i) in page 11 and the hypotheses on cj) and ip, (/>(-, n) is analytic in the deformation 
region and thus, using the inductive hypothesis, the integral equals 

m—k 


^2 / dw k+1 Qp(w k+ i,n k+ i)(f)(w k+ i,n k+1 ) ^ 

n k+1 =1 Z7T1 ^ C o,r-n k u ,k p ,kf>0 ' "' p ' 

k u +k p +kf=k 


k\ 


k u \2 kp k p \k$\ 


x to -vfc +k f [~ dw p dw i V0(w u ,w p ,w i -n u ,n p ,n i ) 

. ™-n fc+1 ( 2vri ) P+ f Jcfxcf (3.16) 


x 'Pq u (w u ] d u )VQ p (w p , uJ { ; n p ,n { )'Pf)(w u , uJ p ,w f ;n u , n p , n f ) 
x II \< II(-irV«+ 1)/2 (l - r n °-)w a p iJj(w p ,w~ p -, n p , n p ). 

a =1 a= 1 

Now we have to deform the Wk+i contour to C^ k+1 . We list the singularities that are 
crossed: 


(a) Since <fr is analytic in Wk+i in the deformation region, the first line of (3.16) only 

contributes the poles of g p (wk+i,nk+i) as specified in (c) in page 11. By (3.1) it 
is clear that computing the associated residues corresponds to grouping (wk+i, rik+i) 
together with the unpaired variables (i.e., turning it into (u;£ +1 , +1 ), keeping (3.8) 

in mind). 

(b) Now we consider the second line of (3.16). The factors involving cj) and ifj are analytic 

in Wk + 1 in the deformation region. As for the factors involving fji, by (d) in page 11, 
the only ones which contribute poles are those of the form f)i(u;)j, Wk+i', rik+i) (see 
in particular (d.ii) ). In view of (3.3), this corresponds to removing (ir^,n^) from the 
group of free variables and adding it to the paired variables together with (wk+i, rik+i), 
which becomes its pair (i.e., the pairs (w^,n { a ) and (u^+i, rifc+i) become (w p ^ +1 , +1 ) 

and (Vfcp + i; n fc p +i)i keeping (3.7) in mind. 

(c) The last two lines of (3.16) do not depend on Wk + 1 , and thus contribute no additional 
poles. 

In addition to the residues in (a) and (b) above, which respectively increase k u and fc p by 
1, we get one more term with Wk+i integrated on the contour Cp k+1 . As in the above cases, 
this can be thought of as grouping (wk+i, rik+i) together with the free variables, this time 
increasing k{ by one. 

The above discussion explains the form taken by the integrand on the right hand side of 
the identity claimed in the proposition (recall that we are passing poles from the outside, 
accounting for additional minus sign in each case). It only remains to verify that the 
inductive argument also yields the combinatorial factor kl/(k u \2 kp k p \k{\) . Consider first 
those terms in the resulting sum which are grouped into sizes (fc u + 1,2A: p , Aif) (which 
satisfy k u + 1 + 2 k p + k{ = k + 1). In terms of our inductive argument, these terms come 
from deforming the last contour when the first k variables have been grouped into sizes 
either (k u , 2 k p , fcf) (so the last one is set as unpaired), (k u + 1, 2(fc p — 1), k{ + 1) (so the last 
one is set as paired, which can be done in k{ + 1 ways), or (k u + 1, 2 k p , kf — 1) (so the last 
one is set as free). By the symmetry of the integrand, and since 


k\ 


+ 


k\ 


k n \2 kp k p \k { \ (k u + l)!2 fc p -1 (fcp - l)!(fc f + 1)! 


(kf + 1 ) + 


k\ 


(fcu + l)!2 fc pA: p !(fc f - 1)! 

(fc + 1)! 


(k n + l)!2 fc pfe p !fc f ! ’ 
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we obtain the desired combinatorial factor on the resulting formula. The other cases (where 
k p or kf are increased by one) can be treated similarly. □ 


Proof of Proposition 3.2. In Proposition 3.3 set <fi = fif 2 and ip 2 = 1)2- By (a) and (b) in 
the list in page 11 as well as the paragraph preceding (3.5), these functions satisfy the 
hypotheses of the proposition. The result follows from the formula given in that result, 
noting that, by (3.2) and (3.4), many of the f 2 factors disappear. □ 


We get now to the crucial point. The only factor depending on x in (3.15) is f 2 (w i a ,n i a ), 
which corresponds to the free variables. As we mentioned, | 1 J 7 J°/' a ' r ~ 2na | < 1 whenever 

f —- f a 

\w a \> t 2n “- Consequently, 

lim hi w ai n a) = 0. 


Now we apply to this to the computation of the limit as x —» oo of p^(f, 2x), which in 
view of Proposition 3.2 amounts to computing the limit of k p ). Since all the sums 

in (3.15) are finite, while all the integrals are over finite contours and the integrands are 
continuous along them, we may take the limit inside and conclude that every term with 
k{ > 1 vanishes in the limit. Using this, and keeping in mind the notation introduced in 
(3.9)-(3.14), we deduce the following 


Proposition 3.4. For every m € Z>o we have 

4%(t) ■= Jim ?££(*, 2s) = ]T - ][ pk , 3£f(fc u , k v ) 

fcu,fc P >o u ' p ' 

k u -j- 2 k p = k 


with 


(?,3»,3r)eA^ >fcp V ; C ° 


w p : n u , n p ) 


ku 


x P 0 u (P u ;n u )p 0 p (h;P;nP)Pf)(P u ,pP;n u ,nP) \w° a 


a= 1 


IP-d 




a =1 


3.3. Simplification of the 1) factors. The next step is to observe that there is an im¬ 
portant simplification in the factors appearing in Vi) (which, we recall, is defined as the 
product of Vi )i and Pfo), due to the emergence of the paired structure. Define 


c(wi,w 2 ;ni,n 2 ) 


(1 — T ni WlW 2 )(l — T n2 WlW 2 ) 
(1 — W\W 2 ){\ — T ni+n2 WlW 2 ) 


Lemma 3.5. For any w\,w 2 € C and n\,n 2 € Z>o for which both sides are well defined, 

i)i(wi,w 2 -n 1 ,n 2 )i) 1 (T~ ni /wi,T~ n2 /w 2 -,ni,n 2 ) 

x i) 2 (T~ ril /wi,w 2 ;ni,n 2 )i) 2 (iui,T~ n2 /w 2 -n 1 ,n 2 ) = T~ niri2 e(w 1 ,iu 2 ;n 1 ,n 2 ). 


Proof. Let hi :2 (wi,w 2 -,ni,n 2 ) = i)i(wi,w 2 \ni,n 2 )i) 2 (T ni /wi, w 2 , n 1 , n 2 ). Noting that the 
first factor in each of the Pochhammer symbols in fq cancel with like factors in i) 2 we get 


hi, 2 (wi,w 2 ;ni,n 2 ) 


[tw\w 2 ] t)qq (r 1+ni+n2 wiw 2 ;T)^ 
(T 1 +ni wiw 2 ; t)^ (t^^wiw^t)^ ' 
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Then the product h\ y2 (wi, w 2 \ n±, n 2 )h\^ 2 ( 1 " Ui /wi,t n2 /w 2 \ ni, 112 ), which is what we are 
trying to compute, equals 


(r l+n ' +n 2 w l w 2 -,T) oo (r 1 ni n2 /(wiw 2 )\t)^ 
(T 1+ni w 1 w 2 -,r) 00 (t 1 ^ 1 / {w\w 2 )\ t)^ (t 1 + 112 w 1 w 2 ;t ) 00 (t 1 ^ 2 /(w 1 w 2 );t ) 00 

Writing w = w\w 2 , the product of the first two ratios above equals 


n\ — 1 

n 


£=0 


1 — T l+1 W 
1 — j yj 


n\—l 

n 


e=o 


1 — t £+1 w 
1 — r ni ~^~ 1 w 




1 — T ni W 


1 — W 


ni — 1 


e=o 


while in a similar way one checks that the product of the last two ratios equals 


n\ — 1 

n 

^=0 


l _ T l-ni-n 2 +£jyj 
1 — q-l-\-£-\-ri 2 yj 


1 - r n2 w 

^ + 7^2 yj 


_ n\ — 1 


n 

1 =0 


-1 


r ni+ri2— i-lyj ' 


Multiplying out the last two expressions gives the result. 


□ 


Now note that, in view of (3.7) and (3.8), r~ n ^/w^ = w%, w a p = T~ n ^/w p and T~ n “/w a p = 
w%. Recalling the definition of the product V\)iP\] 2 through (3.12) and pairing the factors 
coming from the last two products in a suitable way, the lemma implies that 


Vt)(w u ,w p ;n u ,n p ) =Vh(w u ;n u ) J] r-«e«, w p ; n“ n p ) 


x II r (3.17) 

l<a<6<fc p 


Finally, we introduce the change of variables 


w p = t -H z p 


w~P = 


and 


w n = r -H z » 


(3.18) 


(Note that, in view of (3.7) and (3.8), z p and Za P satisfy z p = 1 /Za P , while z “ = cr a ). 
For convenience we also introduce the following notation: given functions g(z , n) and 
h(z\, z 2 ;ni,n 2 ) we write 

g(z,n) = g(T~ 2ri z,n) and h(z 1 , z 2 \ ni, n 2 ) = h{r~^ ni Zl , T~2 n2 z 2 ; ni, n 2 ). (3.19) 

n p 

Under the above change of variables the contour C 0 a is mapped to a circle C 0 t -, 7 of radius 
t~ v centered at 0. Using this and (3.17) in Proposition 3.4, together with (2.20), (2.25) 
and (2.26), and multiplying by the factor ]~[ o r _n “/ 2 introduced by the change of variables, 
we obtain the following: 


Theorem 3.6. For every m € Z>o we have 


E flat( r §mfc(t,0)) = mr ! ^ 
k =0 


flat 

k,m 


(*)= 


E 

Zcu-l-2/cp^/c 


1 !2 fc p fcp! 


3 T(k 


U! ^p)> 


with 
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where 

3m at 0u, k p ) = 


E ! K ^Pf 1 (^,^;n^nP)P0 u (^;n u )P0 p (^;nP) 

Z7T1 pp r kp 
nPlGAP , V ’ JC '0,l 


,n u ,nP)eA^ k ^ 

x II T“ 2n “<c( 2 ;P,zP;nP,nP)e(^ p ,zP;nP,nP) ]^[ r” n >&e(z", z p ; <, n p ) 

l<a<6</c p 


Vb(?-,n«) n \t-i<zZ J](-1)»! t - 5 K) (t-< - 1 )z?UzZ,z?-,<,nl), 


a= 1 


a= 1 


(3.20) 


Co.i is a circle of radius 1 centered at the origin, z a v = 1/za and zf = Ga¬ 


in particular, this proves Theorem 2.5. Note that in (3.20) we have taken the Za contours 
to lie at circles of radius 1 instead of t~ v as specified by the change of variables (3.18). We 
may do this at this point because, although both Q p (za,na) and 0 P {z a p ,n a ) have simple 
poles at Za = ±1, the singularities are cancelled by the double zeros of f) 2 (za, z a p ',na,n%) 
at z v a = l/z% = ±1. 


4. Moment formula: Pfaffian structure 

In this section we will prove Theorem 2.1, starting from the formula given in Theorem 
3.6. The structure of the proof is similar to the one provided in [LDC12] for the case of 
the delta Bose gas. 

In order to unveil the Pfaffian structure lying behind the formula in Theorem 3.6 it will 
be convenient to consider first the case when k u is even. 


4.1. Case k u even. Define the following variables: 

V2a-1 = V2 a = l/z% V2k p +b = <T b , 


(4.1) 


m 2a -i = m 2a = n p , and m 2k +b = n b 


for 1 < a < k p and 1 < b < k u . We will see next how the pairing structure, and in particular 
the simplification presented in Section 3.3, yields formulas which can be naturally turned 
into a Pfaffian. 

A computation shows that 


II ~ e ( z a,zf;nl,nl)e(z a p ,zl,nP,nl) = 

l<a<6<fc p 
and, similarly, 


Tz mb y b - T 2 ma y a 
T ^ ma+ ^ mb y a y h - 1 


(4.2) 


*( z ai z b\Ki n b) = 


A mb y b — t\ m2k P+ a 


U2k p +a 


T2 m “+2 m * k P+ b y a y 2k , b - 1 


(4.3) 


In order to complete the nonlinear Schur Pfaffian (2.32) in J/i,..., y 2 kp +k u we need to 
multiply by some missing factors, yielding 


(L.H.S.(4.2)) x (L.H.S.(4.3)) = 


n 


T 2 m a + 2 m b y a y b - 1 


q- 2 '^'b q /, n~ 2^' a 11 

2kp-\-l<a<b<2k p -\-k u r ^ rz Ua 


k P _im2a-l + im2a,,„ 1 

TT _ 2 y2a-iy2a ~ -l 

, 7 -l m 2a 7 , _ 1 -f^TL2a — 1 7/ . 

a=1 r 2/2 a ~ 2 i/2a—i 


T * mb y b ~ T2 ma y a 

T \{m a +m b )y a y b _ i 


(4.4) 


a,6=1 
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Now consider the factor Vb(z a ',n u ) on the right hand side of (3.20). Since z% = a a , (3.5) 
and Lemma 3.1 yield h(z“, zg; nj), nj)) = (— o- a a b ) n “ ATlb ^ bT? b • Using this and (4.4) 

1—craO-(,rU na+r V 
1 u 1 u 

I "2" ^6 "2 I 

together with the facts that a a = V 2 k p +a, abT i u ~ <TaT i n u = Va<?b sgn(a b n£ - cr a ri£) and 


= 4(r"S-l) ieldg 

r2 m *“ OTo _r3 m i»“-iy 2a _ i (l-(zg) 2 ) 


-(7 6 T^ 


^(2 U ; n u ) x (L.H.S.(4.2)) x (L.H.S.(4.3)) = (-o- a o- fe ) n “ An 6 +1 sgn(<r a n" - a b n b ) 

1 <a<b<k u 


Z%(T n * - 1) 
=1 T5"S( 1 _ ( Z P)2) 


ni 


pf 


T2 mb y b - T2 ma y a 
T \{m a +m b ) yayb _ l 


(4.5) 


- a,6=1 


It remains to recognize the cross-product on the right hand side of (4.5) as a Pfaffian. As 
will be clear in the proof of the next lemma, the factors in that product can be thought of 
as the entries of a degenerate Schur matrix. In the case o a = 1, this identity was discovered 
and checked for small values of k using Mathematica in [LDC12], 

Lemma 4.1. For k € {2,4, 6 ,..positive integers mi,..., mk, and or,... cr*, G {—1,1}, 


(-cr a cr b ) maAmb+1 sgn {a a m a - cr b m b ) = Pf [(~cr a cr b ) ma ‘ Amb sgn(a b m b - cr a m a )] ^ fe=1 . 

l<a<6<fc 


Proof. If any a a m a = CT b m b , both sides vanish, as sgn(0) = 0, so we can assume they are 
all distinct. Switching ( cr a ,m a ) with a different (a b ,m b ) induces a sign change on both 
sides of the identity, hence by making finitely many such switches, we can assume that 


(7; 
form 


= 1, i = 1 

A D 
-D J B 


,..., n, and = —1, i = n + 1,..., k. Our matrix then has the block 
where A is k x k, B is (A; — n) x (k — n) and the k x (n — k) matrix D 


has all entries D ab = 1. Since D is rank one, Pf 


A D 
-D t B 


= Pf [A] Pf [5] by (B.5) below. 


Choose p(h , m) to be a function of h £ 
m is even, and p(h, m) -A 0 as h 
for a, b = 1,..., n and B a ^(h, m) = p( h b’ m b)-p{ h o.,m a ) 
the property that if h a -A 0 in increasing order of m. 


m € Z satisfying p(h, m) -A oo as h -A- 0 if 


0 if m is odd. Let A ajb (h,m) = 


p{hb,mb)+p(h a ,ma,) ^ OT a ’^ " + 


1, 


, k. They have 


(i.e. if m b > m a , then h a 


0 

followed by h b -A 0) for a, b with o a = a b = 1 (i.e. 1 < a, b < n) and in decreasing 
order of m a for cr a = a b = —1 (i.e. n + 1 < a < k) then both A a ^ b (h,m) and B a)b (h,m ) 
converge to (— a a a b ) maA ' mb+1 sgn(cr a m a — a b m b ). By the Schur Pfaffian formula (2.31), we 

have Pf[A(/i,m)] Pf[B(/i, m)] = n l<a<b<n J ^’ a i b(h, m) ]ln+i<a<6<fc B a,b(h, rh). Taking the 
limit in the given order yields the result. □ 


Using this result, (4.5) and the identity t) 2 (za, z a p ; n%, n%) 
(3.20) yields (recalling that a a = y 2 fc p +a and z£ = 1 /z~a) 


(i-Cg) 2 ) 2 

(r»5-(zS)a)(T-"g-(*S)a) 


in 
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CTe{—l,l} feu , 

J2 n a + 2 T, n a=m 


(2vri) 


Sj? n -- 2<< n n 

0,1 l<a<b<kn a<k„.b<kn l<a<b < 


— — T) U T) U 
j- 2 Tl a Tl b 


r «a \ 2 


x 5 P; n u , nP)Pfl u (^; n u )p0 p (^; n?) J] 

< 2=1 

X nf r~ n °(l — (Zg) 2 ) 2 (1-T 

111 J (r "5 - (^) 2 )(r-< - (zP)2) ((z£) 2 - 1) 

2fcp+fc u 

Pf [(-y a yb) maAmb sgn(y b m b - y a m a )]^ 

L J a,6=1 


a=l 


X Pf 


T 2 m t»y fe - T2 m °-y a 

T h^+^)y a y b - 1 


2fcp+fc u 
,l;— 2fc p + l 


(4.6) 

Note that we have removed the restriction that the <r a n“’s be different at this point. We 
may do this thanks to the second Pfaffian on the right hand side. In fact, by Lemma 4.1 
this Pfaffian contains a factor 0 a <6 s S n ( (T fc m 6 — ^a^a), which vanishes if a pair of cr a n“’s 
are equal. Now note that 

5 (Z 7745 (Hz r.\ . T ~ n S l ~ z 2 ) ... = (-^~" /2 -^) no (-^~" /2/ -^) no (P + ^VUr 1+ : V;V L 
0pl ; J0pf /") ) ^ T n_ z 2 ^ T -n_ z 2 'j (- t "/ 2 z;t) ~r"/ 2 /z;r) (tz 2 (r/z 2 

With this in mind we introduce the notation 


u(z,n) = (l-r n )f 1 (r-3 n «;n) 


, \ -In 1 ~ r?1 " 2 

U u (z, n) = t 2 z- 


1 - T n ’ 


(—r n / 2 z-, t)^ (tz 2 ; t)^ 
Up( Zj n) = (-l) n r-3 


-in 1 


: 2 — 1 


(4.7) 


Uap(^,n) = (-l) n T 2 


-In 1 + 2 


z 2 - r 


chosen so that 

fc P 

(2 nd line (4.6)) x (3 rd line (4.6)) = , ^;n p ,n u )Vu u (^-n u ) u p ( 2 >p)t4("») 2 


a= 1 


(4.8) 

(u ap is a certain antisymmetrization of u p , see (4.12), and will be used later on). Of course, 
this decomposition is not unique, but it will turn out to be convenient below. 

Observe that, in view of (4.1), all the y a ' s live in Cop. Therefore we may replace the 
integration in zP over C 0 ^ by an integration in y over C 0 ^ u by introducing suitable delta 
functions. We may similarly sum over mi,... ,m 2 k p +k u € Z>i and get rid of the sum over 
0 i,..., Uk u € {—1,1}. The precise replacement we are going to use, keeping in mind (4.1), 
can be expressed as follows: for a suitable function /, 


E 


1 


(27Ti) fc P Jc*P 
ffe {-lpi^.nPez^^ezJp 0,1 

J2 n a+ 2 J2 n a=m 


dz? f{z p -X p -,n p -,n p -n u ] 5) 


oo If 

^ (27Tl) 2k P+ ku I s-t^kp+ku ^ 11 V2a 

■■■ > ^2fcp+Jfe u = 1 » V a=l 


2kr>-\-k u 


TYl\ ,..., 77 l 2 fcp+fcu 

m H- \-m 2 k p +ku =rn 


_ i-Llm2a-l=m2a 

2/2a 


(^2/a-l + <V + l) 


a=2/c p +a 


X /(?/l, 2/3, • • • ,2/2fcp-i;y2,2/4, • • • ,2/2fc p ; • • • :”l2fcp-i; 

m 2 ,m 4 ,... ,m 2 fc ;m 2 fc + 1 ,..., iR 2 fc p +fc u ; 2/2fc p +i, • ■ ■ ,y2fc p +fc u )- (4.9) 
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The proof of this (essentially algebraic) identity is straightforward. Here, in view of the 
fact that we are working with contour integrals in the complex plane, our delta functions 
implicitly carry a factor of 2vri (which is omitted for notational convenience). 


Remark 4.2. The use of delta functions here will should be regarded essentially as a purely 
notational device which will allow us to express our integral kernels in a more convenient 
way. In most parts of the argument it will not cause any difficulty, with the exception 
of the argument following (4.20) which is treated in pages 26-28. In particular, the delta 
functions will be gone by the time we get to our main result, Theorem 2.1. The reader 
who feels uncomfortable with this treatment can replace 5 Z by a symmetric modifier (p £ (so 
that <p e (x) —> S x as £ —^ 0 in the sense of distributions), and replace the right hand side 
of (4.9) by 


E 


1 


mi ,-,m 2 k p +k_ 

miH-t-m 2 fc p+ fc u =m 


= 1. ( 27Fi ) 


2kp-\-k x 


L 


^2fcp 

° 0,1 


-\-k\ 


dyf (*) 


kp 2,kp-\-ku 

x n^- 1 - l /y 2 °) ^-m,2a-l=m2a n Oftsafoa 

a— 1 a=2kp-\-a 


1 ) + feaiya + 1 ))) 


where /(*) stands for the factor (with the same arguments) appearing on the last two lines 
of (4.9). This is what we do in pages 26-28. 


Using this idea, (4.8), and noting that by the definition (4.1) of the m a ’s in terms of the 
nJJ’s and n%’s 

kp 

n 


r -i (»» 2 


n r- 

l<a<b<kp 


a— 1 

(4.6) becomes 1 

3 W = ;A- 


' n 

a<k u , b<kp 


ii p 
T ~ n >b 


n 

l<a<b<k u 


T 2 n a n b — 


l<a<b<2k p -\-k u 


1 yk 1 f 

2 fcu 2-/ (27ri) 2fc P+ fcu r 2k F 

m 1 ,...,m 2kv+ku =l, v ' 0,1 




dy 


n 


T ~2 r n a m b 


Kp+fcu 
miH-hm 2 fc p+ fc u =m 


l<a<b<2kp-\-k u 


2kp-\-ku 2kp+k u kp 

x Y[ U (y a ,m a ) u u (y a , m a )(6 ya -i + 8 ya +i) JJ u p (?/ a , m a)d y2a _ x - Im2„_i=m2„ 

a=l a=2fc p +l a=l 


X Pf 


T2 m by b - T2 ma y a 
T^+m b )y a y b _ l 


2kp-\-k u 


a, 6=1 


Pf [(-y a y b r aAmb sgn (y b m b - y a m a )] 2 ^ b l 


2kp-\-k u 

2k p +l 


(4.10) 

Now observe that the integration measure ,_ vm u t = m f^< 2k p +ku dy is symmetric 

1 i i 2kp-\-K\x Oq ^ 

under exchanging m a ’s and y a ’s. Let n a € S 2 k p denote the transposition of elements 2a — 1 
and 2a and let T 2 k p be the subgroup of S 2 k p generated by {7r a , a = 1,... ,k p }, which has 
order 2 fc p. Noting every factor in (4.10) is invariant under the action of T 2 k p on the 2/c p -tuple 
((mi, yi),..., {jn 2 k ,y 2 k p )) except for the first of the two Pfafhans, which is antisymmetric 


^There is an apparent singularity in this formula, because u p (y a , n a ) has simple poles at y a = ±1 and y a 
lives in a circle of radius 1. This is resolved by noting that, in view of (2.32), the first Pfaffian on the third 
line has a factor (- r m W 2 y 2a — r m2a_1 / 2 j/ 2 a-i)/(T^ m2Q_1+m2a P 2 i/ 2 a-iJ/ 2 a — 1), which has zeros at these points 
when m 2a -i = rn 2 a and y 2 a = l/f/ 2 a-i- In later formulas it will be less apparent how these singularities are 
canceled (see e.g. Proposition 4.4), but we know that this will be the case by virtue of their equivalence to 
(4.10). 
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under this action, we get 

OO 

3 r(ku,M = * E 


2^11 * (27ri) 2fc p +/Cu /^2fcp+fc u 


dy 


2kp-\-k u 


777-1 v}‘^ , 2fcp + fc u _ 

7771 H-h777 2 /cp + /c u = m 

2/cpH“Ai u 


°o,i 


n 


r ~2 m a m b 


l<a<b<2k p +k u 


X 


u(y a , m a ) u u (y a ,m a )(5y a ^i + 5 ya+1 ) 

a=2/c p +l 


a=l 


(4.11) 


2/cpH - A;u 


X Pf [{-yayb) maAmb sgn (y b m b - y a m a )] a ^ = . 


2fc p +l 


Pf 


T 2 m »y b - T2 m “y a 

r | (m a +m h ) ya y b _ 1 


- a,6=1 


X] S S n ( 7r ) II^C^C 2 ' 1 - 1 )’" 1 ^ 2 ' 1 - 1 ))' 5 : 




a= 1 


y-n(2a — 1) 


r(2a) 


^7r(2a —1) ^7r(2a) * 


Observe that the product of the delta and indicator functions in the last line is invariant 
under exchanging any pair of indices of the form (2 a — 1, 2a), a = 1,..., k p . On the other 
hand, we have 

i m l + y 2 


u P (y, m) — u p (l/y, m) = (—1)”V 2 ? 


1 - y“ 


= Uap (y,m), 


(4.12) 


where u ap was defined in (4.7), and thus it is not hard to see that 

1 fc P 

(last line (4.11)) = -j- JJ U ap (l/2o-l, ™2a-\)&y aa - l -\/ y2a lm aa -i 


-l=m 2a 


0=1 


Using this and repeating the procedure we just used of using the symmetry of the integra¬ 
tion measure to sum over permutations, the last formula for 3m (& u , fc p ) can be rewritten 
using the replacement 

kp 

(last line (4.11)) X 2^2*0! II U ap(2/,r(2a-l), rn a ( 2 a-l))S y - 1 1 r 


2^p(2fcp)! 

a£S 2kp a=1 


V(r{2a — 1) 


r i 2 k 

Since, by (4.12), the matrix [u a ,p(ya,yb\rn a ,m b )l ma=mb 6 ya _ 1 / yb l c tf b \ a ^ 1 is skew-symmetric, 
this sum equals pk’y. times the Pfafhan of this matrix. Thus, using in addition the identity 
Ha=i A a Pf[Al(a, b)] 2 ™ b=l = Pf[A a AfeA(a, b)]^ l b=1 (see (B.4)), the last sum can be written as 


t (2 a) 


' l a(2a-l)— m a(2a) ' 


1 k ' 

1 V Pf 


2 2fc r (2fep)! 


1 2 fc D 


4u ap (y a , rna)lm a =m b fiy a _L la^fe 


Vb 


- a,6=1 


(the introduction of the factor 2 2kp in front of the Pfafhan is for later convenience), and 
now using the same Pfafhan identity again we deduce the following 


Proposition 4.3. Assume that k u is even. Then 


3r(*u,*p) = 


v 


E 


1 


L 


2fcu+2fcp/'2T ji f27rij 2fc p+ fc u /„2fc p +fc u 

V p ; mi,...,m 2 fc p+ fc u =l, V ’ J °0.1 

mil-t-m 2 fc p+ fc u =m 


dy 


n 


T ~2 m a mb 


l<a<b<2kp-\-k u 


2/cp+fcu 


u(y a ,?Tla)Pf 


a=l 


- 2kp-\-k u r\i 

T2 mb Vu — T2 ma ll 

Pf 4u ap (y a ,m 0 )l ma=mb 8 il a ^ b 


T 2^na+m b ) ya y b _ 1 


a,6=1 


Vb 


a,6=1 


X Pf [(-yayb) maAmb sgn {y b m b - y a m a )u u (y a , m a )u u (y b , m b ) 


2kp-\-k u 


x(6y.- 1 +6,„ +l )(6 tb -l + Sy b+ l)YX 


2kp-\-l 


(4.13) 
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Note that the product of delta functions in the last line of (4.13) poses no difficulty when 
0^6 (since they involve different variables), while for a = b the prefactor sgn {y b m b — y a m a ) 
vanishes and thus the whole factor is interpreted as 0 . 

Recalling now the definition of pju^(i) in Theorem 3.6 we have 


flat 

k.m 


(*)= 


E 


fcn ,fcp>0, fcuH - 2/cp— k 


2 k (2k p )\k v 


r(*)> 


(4.14) 


where (*) stands for the expression starting from the first sum in (4.13). Note that in this 
expression the only factor which depends directly on k u or k p (as opposed to k u + 2k p = k) 
is the second Pfaffian in (4.13). Remarkably, this Pfaffian can be combined with the third 
Pfaffian in (4.13) and with the sum over k u and k p in (4.14), thus yielding an expression 
which only depends on k (since, furthermore, the combinatorial prefactors also combine 
into something which only depends on k). This is done using Lemma B.2, and the result 
is the following: for m > k > 0 , k even, 


4t(t) = 


{~l)z k 
2 k kl 


E E 


771 lH-1-771/2=777 


(27ri 


;\k 


r. k 

I dy n T-^^Wuiya 

Jc 0 ,i l<a<b<2k p +k u a=l 


m„ 


x Pf 


(- 2 /a?/ 6) maAm6 sgn {y b m b - y a m a )(S ya -1 + S ya+ i)(S yb - 1 + S yb+ i)u u (y a , m a )u u (y b , m b ) 

1 k 


T 4u ap (?/ a , n 7 . a )l ma=r7l(j 5 _ l_ 1 a^b 


Vb 


- a,b= 1 


Pf 


A^ya - A m »y b 


T - 2 {m a+ m b ) yayb _ 1 


- a,6=1 

(4.15) 


Notice that in the last Pfaffian we have flipped the sign of the argument, which yields the 
additional factor (— l ) fc / 2 in front. This will be convenient later. 


4.2. Extension to general k u . To handle the case when k u is odd we use a standard 
trick (which goes back at least to [DeB55]) to extend Pfaffians to matrices of odd size (the 
same idea is used in [LDC12]). Consider the dummy variables 


U 2 k p +k u +i — C and rri 2 kp+kp +1 — 0 . 
A simple computation (together with (2.32)) shows that 


n 

l<a<b<2kp-\-k u 


Vb - Va 

VaVb ~ 1 


n 

1 < a< b< 2k p + k u +1 


Vb Va 
VaVb ~ 1 


Pf 


Vb - Va 

_VaVb - 1. 


2fc p +/c u +l 

a,6=l 


and, similarly (using Lemma 4.1), that 


JJ (-cr a (Tb) maArnb+1 sgn (a a m a - cr b m b ) 

2kp+l<a<b<2kp-\-k u 

= Pi[(-cr a cr b ) maAmb sgn (a b m b - a a m a )\ „^ 2 fcp+i' 

These formulas can be used directly as replacements for (4.5) and Lemma 4.1, and lead to 
a formula which reads exactly like (4.10), except that the matrices appearing inside the 
two Pfaffians are augmented in this fashion. This can then be extended directly to (4.13) 
and further to a version of (4.15). We leave the details to the reader and simply record the 
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result: for odd k u we have 


/_"I 5 §(fc+l) 

k ,m{ ) 2 k kl 


mi 

mi H- \-m^=m 


l<a<b<2kp+k u 


o° 1 /■ 

e e pi/,« n 

k 

x n u(y a ,m a )Pf 


a=l 


A U 
-U J 0 


B 1 
-1 T 0 


(4.16) 


where A and B are, respectively, the matrices appearing inside the first and second Pfaf- 
fians in (4.15), U is the vector with entries XJ a = u u (—1, m a )6 ya+ i — u u (l, m a )5 ya _i and 1 
represents a vector of ones of size k u (which is formally obtained from (4.15) by augmenting 
the matrices in the Pfaffians using the dummy variables and computing their last row and 
column). 

Now Lemma B.l with V = 1 together with the identity (— l) 2 fc ( fc+1 ) = (—1) 2fc lfce2Z + 
(— l)z( k+l) lk & 2 Z+i allow us to write a version of (4.15) and (4.16) which is valid for k even 
and odd: 


Proposition 4.4. For any 0 < k < m we have 

f_l)i fe ( fe + 1 ) EE. 


r_ “ /■ i 

A k %(t) = { —L ^— L<® n >(».,».>«[*(»;«)] 

^0,1 l<a<6<2/c p +fc u a=l 


mi 

mi H- 1 -m^=m 


^ ^ us + k u 

where was defined in Proposition 3-4 and 

1 k 


_ [Ki,i(y a ,yb',m a ,m b )] , ± [K h2 (ya,yb;m a ,m b )] b=1 
Jl\ \l)\ TTlj — jk p 

[KiAVbi Va ; rn b , m a )\ a b=l [K 2 , 2 {ya, Vb\ m a , m b )\ a b=l _ 


with 2 


-^1,1 idJai Vb'i TBai m b) — 4u a p(t/ a , 771 a )l ma=m(j (5_ L^-a^b 

u Vb 

+ (-y a yb) maAmb sgn (y b m b - y a m a )u u (y a ,m a )u u (y b ,m b )(Sy a ^ 1 + d Va+ i)(5 Vb - 1 + <5j, i+ i), 
Ki, 2 (ya,yb',m a ,m b ) = u u (-l ,m a )S ya+ i - u u (l, m a )5 ya _i, 


K 2,2 {y a ,yb-,rn a ,m b ) = 


y a - T 2 °y b 

T ±(m a +m b ) ya y b _± 


(4.17) 


where u, u ap , and u u are given in (4.7). 


This result together with (2.20), (2.25) and Proposition 3.4 yield the following explicit 
formula for the exponential moments of flat ASEP: for any m € Z>o we have, with K as 
in Proposition 4.4, 


E : 


flat 


r ±mh(t, 0) 


= m. 


■E 

k =0 


(_l)ifc(fc+!) 


dy 


2 k k\ J C i 

= 0,1 

raiH- \-mi £ =m 

k 

r -|mam 6 JJ u (y a) ma ) Pf[A"(y; rh)\ . (4.18) 


X 


l<a<b<k 


a= 1 


o 

In order to avoid an even heavier notation, we are slightly abusing it in the definition of 
Ki'i(ya,yb',m a ,mb), which depends on a and b explicitly in addition to y a , yb, m a and mb- This will 
become irrelevant in later version of our formulas, see footnote 4 in page 29. 
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At this point we can get rid of the delta functions inside K , by rewriting (4.18) in such a 
way that we can take the y a integrations inside the Pfaffian. To this end, let p a = , 

so that )/2 ^ . Note that 9i(p a ) > 0. The key will be to use the identity 

Pa ~ Pb ' = [ d\ a d\ b e- x “P“~ x >P» sgn(A b - A a ). (4.19) 

PaPbVPa + Pb) J{ R> 0 )2 

In order to use it, we first use the Pfaffian identity (B.4) to rewrite (4.18) as 


E h 


T \mh(t, 0) 


= m- 


HO -j VJ p 

! E a e n r-i—n-o*..™.)*. 

k=0 v ' °o,i l<a<6</c a 

miH-hw,fc— m 


X (_l)|fe(fe+ 1 ) Pf 


[-^ 1,1 (y«, m a , m b )\ k a b=1 [±K lt2 (y a , Vb] rn a , m b )\ k b=1 


r [^K\ 2 {y b ,y a -m b ,m a )] ab=1 


Pa-Pb 


PaPb(j>a+Pb)l a,6=1 


Now using (4.19) and Lemma B.3, which is a certain Pfaffian version of the Andreief (or 
generalized Cauchy-Binet) identity [And83] (more precisely, we use a version of the lemma 
with the matrices inside the two Pfafhans transposed), with respect to the integration in 
the A a ’s, with 4> a (X a ) = e ~ XaPa , we get 


E : 


flat 




m ^ oo If f 

^ (2^ij* L^L), 

k =0 mi,...,m k =l y °o,i 'P K >ol 


= m T ! > — 

miH-h m k =m 

->*aPa | | t-— 3 rn a mb 


d\ 


J u(y a ,m a )p, 


a= 1 


n r- 

l<a<6<fc 


(4.20) 




I" [Ki,i(y a ,yb;m a ,m b )] k ab=1 [Ki t2 (y a ,y b -,m a ,m b )] k ab=1 

T 1 T -| r 1 fc ? 

_-[-^i, 2 ( 2 / 6 ,ya;"r 6 ,^a)J a>6=1 [ sgn(A b - A a )J a b=l 


We would like next to interchange the y a and X a integrals. The difficulty we face is that 
the y a contours pass through singularities of the factors u ap (y a ,m a ) in the integrand, and 
after interchanging the integrals the zeros of the second Pfaffian in (4.11) (see footnote 1 in 
page 22) are not there any more cancel them. This means in particular that the y a integrals 
that result after the interchange will have to be regarded as principal value integrals. In 
order to justify the interchange of limits let us restrict ourselves to the case k = 2. The 
general case can be justified similarly but the notation becomes much heavier. 

When k = 2 we may use the Pfaffian formula 


Pf 


0 a b c 
—a Ode 
-b -d 0 / 
-c -e -f 0 


to write the double integral in (4.20) as 


af — be + dc 


-1 

4(27ri) 2 



dXu(yi,mi)u(y 2 ,m 2 )pip 2 e Xipi X2P2 r z” 11 ™ 2 


x 


dCi ! i(yi,y 2 ;m 1 ,m 2 )sgn(X 2 


Ai) - dC l2 (yi, yp, m 1; m 1 )K h2 (y 2 , y 2 ; m 2 , m 2 ) 


+ Ki ,2 (yi,y 2 ',mi,m 2 )K 1>2 (y 2 ,y 1 -,m 2 ,mi) 


The last two terms in the bracket pose no difficulty. In fact, recalling that K\ t2 (y, y'\ m, m') 
only depends on y and m, each of the two terms involves products of delta functions 
depending on different variables, and one can check directly that the yi,y 2 integration 
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(which corresponds to summing the values at ±1) can be done before or after the Ai,A 2 
integration. 

Now consider the first term in the bracket above and call I\ the associated integral. 
Integrating Ai,A 2 first and using (4.19) yields 

h = A , 0 [ dyu(y 1 ,m 1 )u(y2,m 2 )T~^ mim2 ^—^R\ 1 (y 1 ,y2-,rni,m2). 

4(2tti r Jc$ tl P 1 +P 2 


K\ \ is defined in (4.17) and has two terms. Call I\y and I\, 2 the integrals associated to 
each (so that I\ = I ip + Up). As before, the second integral poses no difficulty for the 
interchange of X a and y a integrals, so let us focus on Up. Here, thanks to the prefactor 
P ^~ P2 , the integrand has no singularities on the y a contours and thus the integration simply 
yields 


r "I f , , w i n / . ,i 

h,i = — / dyu(y,mi)u(-,mi)u ap (y,mi)uap(-,mi)—- ——-j -- 

2ti Jco,i y y i) 

with p(y,m) = \^Sr y - 


T 




mi=m2 


(4.21) 


Next we let let Up be the analogous integral where we now integrate first in y. In this 
case the integrand has singularities along the y contour, so we will have to regard this as 
a principal value integral while dealing with the delta functions at the same time. More 
precisely, letting Cf } , be portion of the contour Cop which lies outside of the circles of 
radius 5 centred at ±1, and letting p £ be a suitable modifier, we define Up as 


In = lim lim 


-1 


e <5i, a 2 —>o (27ri) 2 7 (r>0 ) 


[ dX [ 

R>n) 2 JCi\ 


0,1 X ( A ),1 


dy u(y 1 ,mi)u(y 2 ,mi)u ap (yi,mi)u ap (y 2 ,mi) 


xr 2 m ipip 2 e Aipi A2P2 sgn(A 2 - Ai)^ e (?/i - l/y 2 )- (4.22) 


Using Cauchy’s Theorem, each of the 5 a —> 0 limits equals the integral over a circle Cop 
of radius slightly smaller than 1 plus half the residues of the integrand at y a = ±1. Using 
this to compute the two limits (one after the other) leads, after some calculation, to 


Up = 


— lim , 
e-s-0 (27ri 


w. 




d\ dy t 2mi u(y 1 ,rn 1 )u(y2,rn 1 )u ap (y 1 ,rni)u ap (y 2 ,rni) 

(M>o) 2 J (C 0 ,i) 2 

x pip 2 e~ XlPl ~ X2P2 sgn(A 2 - Xfippfiyi - l/y 2 ) 


— lim 

£— 


V — / dx f 


.1^2 


cre{-l,l} 


dX I dy t 2 m i u (y,mi)u(fj, mi)u ap (y,r?7.i)o-2(-l) 
(K > 0 ) 2 JC 0 ,i 


m i 


x T~5 mi p(y, m 1 )p(a 2 ,m 1 )e- x ^y^- X2P ^ m ^ sgn(A 2 - X\)ip e {y - o) 

— lim- V [ dXT~ 2 m iu(ai,mi)u(ao,mi)cria 2 T~ mi 
£ ^° 4 w ,u%o) 2 

x p(cri,mi)p(fJ 2 ,mi)e _Alp(,71 ’ mi) ” A2p(<72 ’ mi) sgn(A 2 - Ai)<£ £ (op - cr 2 ). 


At this point we can perform the X a integrals using (4.19) (note that in the first two terms 
there is now no difficulty in interchanging the X a and y a integrations). Starting with the 
last one, it yields a factor ] i when op = op this factor vanishes, while when 

<7 1 o 2 the limit linp^o^(op — op) vanishes. This means that the third limit above is zero. 
One can check similarly that the second limit vanishes, by performing the A a integrations 
to obtain a factor an d th en computing the e —>• 0 limit to yield evaluation 

at y = o, for which this factor vanishes. Performing the X a integrals now in the first term 
and taking e —>• 0 shows that Up equals the right hand side of (4.21) with the contour Cop 
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on the right hand side replaced by Cop. But since the integrand has no singularities along 
Cop, Cauchy’s Theorem implies that Zip = Zip as desired. 

The conclusion of all this (and its extension to all k > 1) is that we may interchange the 
A a and y a integrations in (4.20), provided we mollify the delta functions and we interpret 
the y a integrals as principal value integrals as done in (4.22). Now observe that in the 
Pfaffian appearing in (4.20), the (1,2) entry only depends on y a , the (2,1) entry only 
depends on y b , and the (2,2) entry depends on neither of the two. Hence we can use 
Lemma B.3 again, this time to bring the y a integrals inside the Pfaffian. Let us first define 


o(A ai y a ,m a ) = p a e XaPa u{y a ,m a ) 


il _ T m a /2 -a a '^ +t 


— —(1 — T yria T m a l2 1 r ^ a c l+r m a/ 2 


l+r 


-m a /2 


Va 


1 + r m “/ 2 y a 

( T~ mn/2 ?/; t)^ (r 1+m ^ 2 ;T) c 
(_ r m./2 ya . r j ('Ey 2 ; t)^ 


Wo l+ T m “/ 2 Wo 


(4.23) 


where we have used the definition of u in (4.7) (and of jp in (2.23)), and observe that the 
middle line of (4.20) can be written as r -™ 2 / 4 n„T 3ma,:, (A a , Va, Tn a ) (recall that Yla m a = 
m). Define also 3 


Kff(X a ,X b ;m a ,m b ) = 1 „ / dydy'T^ m ^ +m b^{X a ,y,m a )v(X b ,y',m b ) 

2 ( 2 vri) Jc % 1 

x Ki,i(y,y';m a ,m b ) 

If 12 

= 1 ma =m b —f dyr2 ,n “t)(A a ,y,m a )t)(A fc , l/y,m b )u ap (y,m a ) 

70 Jc 0 ,i 

+ ^ E (-^) m “ A “ i+1 sgn(^-™ 0 ) 

o-.o-'ef-ip} 

x r l( m 2+ m 2)-|(m a+mi) ) D ( Aaj ^ m a )d(Afe, a', m b ), 

Kff(X a ,X b ;m a ,m b ) = ^-r [ dy r* m “t)(A a , y, m a )K^ 2 {y,m a , -) 

4m JCo,i 

= -\ E r>'-^0(A a ,a,m a ), 

cre{-ip} 

K 2 ^(^a, Afe; m a , m b ) = - sgn(A fe - A a ) 

(4.24) 


where u ap is given in (4.7), 0 is given in (4.23) (the —’s in denote the fact that this 
kernel does not depend on those arguments), and we have used the fact that for a € {—1,1} 
the function u u (cr, m) (defined in (4.7)) satisfies u u (cr, m) = r ~~ Z = crr~~. Then, 
by virtue of Lemma B.3 again, we get 


E : 


flat 


-^mh(t, 0) 


_ i™2 

= m r \T 4 


1 


E 


k =0 ^ ^(®>o)' 

miH-h m^=m 


dX 


x Pf 


[it'| i ! a 1 t (Aa,A( ) ;m a ,m 6 )] ^ ft=l [A'f i | t (Aa,A(,;mo,m 6 )] ^ 6=1 
- Aa;mf,,m a )] a i=1 [R'f3 t ( A a>A{,;m a ,m i ,)] ^ fc=1 


(4.25) 


(Note that we have used the fact that Pf[aH] = a n Pf[H] for a skew-symmetric 2 n x 2n 
matrix H and a scalar a. See (B.3).) 


^The notation f in the second line of this formula indicates that this is a Cauchy principal value integral. 
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The formulas for and follow directly 4 from the definitions of d, K\,\ and 
with the nuisance that we have to justify writing a principal value integral in Kf*f. This 
involves writing out explicitly the result of the interchange of the A a and y a integrations 
in (4.20) using modifiers and principal value integrals in a similar way as (4.22), then 
moving the y a integrals inside the Pfaffian, and finally showing that the limits involving 
the mollification and the principal value integrals can be taken back inside the Pfaffian 
to yield the kernel The argument is very similar to the one we used to justify the 

interchange of limits itself, so we omit it. 

The factor (—l) 2 fc ( fc +i) ; n f ron t of the Pfaffian in (4.25) now allows us to write the 
formula in the form of a Fredholm Pfaffian. In fact, let AT flat be the 2x2 matrix kernel 


iP flat (A a , X b ;m a ,m b ) 


K^(X a ,X b ;m a ,m b ) 

-Kfff(X b ,X a ;m b ,m a ) 


Kf^(X a ,X b -m a ,m b ) 

K^{X a ,X b ;m a ,m b ) 


(4.26) 


The matrix ^ flat (A a , X b ] m a , m b )\ ^ b _ l differs from the one inside the above Pfaffian by a 
permutation of the rows and columns which has sign (—l) fe ( fc—1 )/ 2 ) and thus we get 

“ ("!)* 


E : 


flat 


T \mh(tfi) 


= mJ.T~> 2 


£ 

k =0 


k\ 


V / dX Pf 

=1 J(M>o) k 


K &at (X a ,X b -,m a ,m b ) 


i k 


mi~\ - \-rrik=m 


. a,6=1 

(4.27) 


This gives Theorem 2.3. 

As a corollary we get a bound on the moments of r^'^ 0 ), which will be useful later on 
when we form a generating function for 


Corollary 4.5. There is a c> 0 such that for all rn € Z>i, 


E 


flat 


T \mh(t,0) 


< cm T W 


-im 2 


(4.28) 


We note that for fixed r we have m T \ < c(l — t) m (this follows from the q -analogue 
of Stirling’s formula for the q-Gamma function), and hence the above moment bound 
is actually of the form cr"4 m - 0 (m). B e f ore getting started with the proof we need to 
establish an estimate on certain ratios of g-Pochhammer symbols. 


Lemma 4.6. For any r/ > | and r £ (0,1) we have 

3 (-t-" /2 2 ;t) 00 (t 1 +”zV), 


sup T 
|*|=1 


T/n 


{—T n / 2 Z]T^ oo (tz 2 ; t) c 


< oo. 


Proof. By definition of the g-Pochhamer symbol, 

(-r-"/ 2 *^)^ {P+ff£yr)^ rrn-1 1+P-*/£ 

(-r"/ 2 * 2 ;r) (tz 2 ;t) IL =0 l-r ^ 1 * 2 

It is not hard to see that the absolute value of each factor is maximized for 2 in the unit 
circle at z = 1, whence 


(- T ~ n / 2 z-,t) (t 1+ti z 2 ;t) 
(-W 2 2 2 ;t) (tz 2 \t) x 

v / OO 


/ mToP+r e ~ n/2 ) 


<^n,="o(i+^ /2 )- 


Expanding the last product and bounding each term by the biggest one, which is given by 
n[=o 2J T l ~ n / 2 ~ t - ” 2 / 8 results in the bound 


(_ t -"/2 z-t) (t 1+ "z 2 ;t) 

_v_' 7 00 _V_' / 00 

(-W 2 22;t) oo (tzV)^ 




c 


2 n T ~ n2 /^ 


4 Observe that we have dropped the factor l a jtb in the second term of the second equality in the definition 

of Ki f. We may do this because this term vanishes when A a = At, and m a = rrn,, as can be checked using 

the antisymmetry of (to, m) 1 —> sgn(m — m). 
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for some C > 0, which implies the result. 


□ 


y 


Proof of Corollary f.5. Recall that u(A a , y, m a ) contains a factor e~ XaPa withp a = 
and note that $ft(p a ) > c\ for some ci > 0, uniformly in y € Co,i and m a € Z>i. Having 
chosen ci in this way, we use (B.4) to factor out e from the Pfafhan: 

l 


Pf 


AT flat (A a , X b -,m a ,m b ) 


it 


- a,6=1 


n 

a= 1 


e ~Cl\a T ± m l pf 


A' flat (A a , A b ;m a ,m b ) 


- a, 6=1 


where 

JP flat/ (A a , A&; m a ,mb) = 


e c 1 (Aa+A ( ,) T -("*a+"‘6)/ 16 JT lil (Aa,A 6 ;m a ,m ( ,) e c l Xa T~™ 2 a/ le K li2 (\aM;m a ,m b ) 


- e c i x b T m 6/ 16 A'i i 2(A6,A a ;m(,,ma) K2,2(^a,\ b ;m a ,m b ) 

The point of rewriting Lf flat in this manner is that the entries of iL flat (Ai, A 2 ; mi, m 2 ) 
are uniformly bounded for Ai,A 2 > 0 and mi, m 2 € Z>i, say by a constant C 2 > 0. 
This can be checked straightforwardly for most of the factors involved, except for two 
minor issues to address. First, d(A a ,y,m a ) contains a factor of the form g(y,m a ) '■= 

(- T-^-aPy-r ) (r 1+m “y 2 ;r) 

t T m a /2 yT \ (tj/ 2 ;t) 

\ / OO 


L , which is not bounded in m a . Nevertheless, t>(A a , y, m a ) appears 


next to r m “/ 4 , which together with the factor r - " 1 ®/ 16 introduced above gives a factor 
r 3m “/ 16 . Thanks to Lemma 4.6 this factor is enough to balance the growth of g(y, m a ). The 
other term which is not seen directly to be bounded is the one coming from K\ 1 involving 
the principal value integral l ma=mb dy t)(A 0 , y, m a )t>(\ b , 1/y, m 0 )u ap (j/, m a ). But in 

this case the uniform bound follows from the fact that the integrand can be expanded 
around ±1 as jf- + 0(1). 

Using the fact that iL flat is uniformly bounded, the identity Pf[A] 2 = det[A] and 
Hadamard’s bound we deduce that |Pf [/L flat (A a , Af m a , mb)] ^ ft= , | < k k / 2 c k 2 n„ e -ClAo T m “/ 16 , 


and thus 


' (®>o)' 


d\ Pi[K^\\ a ,\ b -m a ,m b )] k ab=1 


< 


(C 2 /C\) k k k / 2 JJ T !6 


(4.29) 


Now we use this bound in (4.27) to deduce that 


E 1 


flat 


T |mh(t,0) 


OO 


<m T \r X J2 ^ Z. 

k=0 m 1 ,...,171^=1, a 
miH-1- m^—m 


n 


T i6 m a c k k k / 2 < dm T \r P 


for some c, d > 0 which are independent of m (using Stirling’s formula). 


□ 


5. Generating function 

We are in position now to turn this into a formula for a certain generating function 
of flat ASEP. We will consider the following generalized ^-exponential function: for fixed 
(€C with |£| < 1, 

00 1 

ex Pg (z;£) = ^— 
k =0 Kq ' 

Note that for |£| < 1 and x £ C this is function is analytic in both variables. On the other 
hand, for fixed £ with |£| = 1, this function is analytic in x for \x\ < 1. Choosing £ to be 1, 
g 1 / 2 and g 1 / 4 yields, respectively, e q (x), E q (x) and exp g (x) (although note that, for £ = 1 
it only yields e q (x) restricted to \x\ < 1), see (2.5) and (2.6). 

The next result gives a formula for the expectation of exp T (£r 2 /l ^ ,0 A £) for |£| < r 1 / 4 . 
The restriction on |£| comes from the fact that we will use our moment bound (4.28). 
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Theorem 5.1. Let £, £ € C and assume that either |£| < r 1//4 or |£| = r 1 / 4 and |£| < r 1 / 4 . 
Then 


E 


flat 


'exp T (Cr3Mt.o) ;e) l = £(-l)' 


/c=0 


1 

fc! 


mi 


1 f 

,5*=i ( 2iri ) k J<*> 0 ) k 


d\ 


1 k 


. a,b= 1 


(5.1) 


x ^E a ^ r -i(E a m a )^(E 0 m Q ) 2 -E Q m apf AT flat (A a ,A fe ;m a ,m 6 ) 

with A' flat as in (4.24). Moreover, the series on the right hand side is absolutely convergent. 
Proof. Using the definition of exp T (-;£) and interchanging sum and expectation we have 


E : 


flat 


exp T (£r^ (t ’ 0) ;£)] = 


k =0 


T \kh{t, 0) 


(5.2) 


Here, in order to justify the use of Fubini’s Theorem, we are using Corollary 4.5, which 


gives 


jgflat 


T 2 


0 ) 


< ck T W ~ k2 / 4 , and implies that the double sum on the right hand side 
is absolutely summable under the stated conditions on £ and £. 

Now rewrite the right hand side of (4.27) as 


oo , . e oo 

E 

£=0 mi.m<=l, 


r 4(E a m a ) 2 


<2(m) 


miH-|-m,£=fc 

with Q(rh ) = f^ >0 y d\ Pf [A" flat (A a , A&; m a , m&)] ( f t 6=r To get a formula for the right hand 

side of (5.2) we multiply this by -^C k f k ^ k ~ l \ sum over k = 0,..., N and then take N —>• oo 
to get 

/ I \£ c*- 1 

l ' _i ' £Ea m « r -l(Ea m “) 2 £(Ea m “) 2 -Ea m «Q( ? A). 


E 


flat 


exp r (CT2 h(t ’ 0) ;£)l = lim Y 

AT—inn L -* 


N^OO Z J t\ 

£=0 mi,...,mg=l, 

m\-\ - \-m£<N 


We claim that the limit can be taken inside the first sum, which yields 

OO / 1 \£ OO 

flat //-^bh(t,0 ). \ A v 


E : 


exp T (Cr^( t ’°);£)j = X] C E “ m “^3 ( E B ™a )2 ,t ( Ea^ )2 -E Q m a g(^) 


£=0 


and proves the result. To see that the limit can be taken inside it is enough to show that 
the right hand side above is absolutely summable. But this follows from (4.29), which gives 
\Q(m)\ < c¥' ,2 r^« m “/ 16 . In fact, using this bound we get 


OO ^ OO 


Y~ |C|Ea^ r -|(Ea^) 2 |£|(E a "la) 2 -Ea^|Q( ? A)| 


t\ 

=0 = l 


oo 

sE 


c(/ ^ - , E a m ar -i(E a m a ) 2 + 1 ^E a m 2 |y|(Ea m «) 2 -Ea 


£! 


E ici ] 


; ier 


£=0 " mi,...,me=l 

which is finite, as desired, by Stirling’s formula and our assumptions on £ and £. □ 

We focus now on the case £ = r 1 / 4 of Theorem 5.1, which yields the exp T -Laplace 
transform. In this case the factor £(Ea m U j n f ron ^ Q f the Pfaffian in (5.1) cancels exactly 
with / 4 ; an d hence the sums in m a can be brought inside the Pfaffian. To this 

end, and in view of the definition of A" flat in (4.24), we define, for £ € C, 

Aff c (A a , Aft) AfJ’ C (A a ,A fe )‘ 

A"i fl |’ c (A b ,A a ) <2 t,C (A a ,A ft ) 


K^(X a ,X b ) = 























A PFAFFIAN REPRESENTATION FOR FLAT ASEP 


32 


with 


A'ff’ c (A a , X b ) = V — f dyT^ m 2 C 2 m v(K,y,m)v(\ bl l/y,m)u ap {y,m) 

El mJ Co, 

oo 

+ 2 E E (-cra') mAm ' +1 sgn(a'm - am) 

m,m'= 1 cr,cr'G{ —1,1} 

X T I(m 2 +(m') 2 )-i(m+m')^m+m' t ,( Aa; ffj m ) D (A fe , a', m'), 

1 OO 

^A) = -jE E r^ 2 -H m t)(A a ,a,m), 

m=l itG{—1,1} 


< 2 ’ ? (Aa, Afe) = isgn(A b - A a ), 


(5.3) 


where u ap , and u u are given in (4.7) and 0 is given in (4.23). Then using Theorem 5.1 and 
Lemma B.3 we deduce the following 


Theorem 5.2. For |£| < r 1 / 4 , 


E 


exp T (Cr§^ 0) )] = ]T 
fc >0 


k\ 


( R >o)' 


d\(-l) k Pf[K^(X a ,X b )] k ab=l 


k 

. a,6=1 


= pf 


J - A' flat ’ C 


- i 2 ([0,oo)) 


Finally we are ready to prove Theorem 2.1, for which we need to check that the kernel 
given in (5.3) coincides with the one given in (2.12). The proof amounts to what can be 
considered as an ASEP version of the “Airy trick” commonly used in the physics literature. 
In fact, the representation of t>(A, y, m) through an (inverse) double-sided Laplace transform 
(see (5.7) below) is an exact analog of the identity 


3 /s _ 


f 


dw Ai (w)e u 


(5.4) 


for m > 0 used in [LDC12] (see (42) and the identity right before (137) in that paper). 
A crucial difference, though, is that in [LDC12] (as in other situations were the “trick” 
is applied), the authors are dealing with a divergent series, and the Airy trick is used to 
resum it through an illegal interchange of sum and integration, whereas in our case such 
an interchange will be fully justified. 


Proof of Theorem 2.1. For simplicity we will write K instead of iP flat A throughout this 
proof. The only two pieces of the kernel for which we need to show that the two definitions 
coincide are K\ t ± and K 12 , and we will begin with the latter. Throughout the rest of the 
proof we fix /? > 0 as in the Section 2.2. 

For C ^ E>o define 


A,y;C) 


(-cr 


-~hs 2 -b 


(n/ 2 ; r ) c 


(1 -r) 


(r l+s y 2 ] t) c 


-f(A ,y,s) 


(5.5) 


Recalling the definition of D in (4.23), this definition coincides with the one given in (2.8). 

(t't) 

Moreover, for m € Z>i and a € {—1,1}, since (i„ r )m(’ T i+m. r ) — = m r !■ if satisfies 


12 1 

sj-jm — -^m /-m 


( — — — 771 

i™C"MA ,a,m) = (-l) m T * 4 ^(m,A,u;C)- 


rn- 
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Therefore 


1 (^ — -\-^^TTl 2 — — 771 

K h2 (X a ,X b ) = --Y: ^2 (~ l ) m ~ - —- V>("b Xa,cr-,()- (5.6) 


771,-! 


m= 1 cre{-l,l} 

Now we introduce the inverse double-sided Laplace transform of ip (in s), given by 
V>(w, A, y; C) = 77 -T f ds e SUJ ip{s, A, y; C). 

^ 7rl JiR 

This integral converges thanks to the factor r~d s2 which we introduced in ip, and we have 
the inversion formula (see [Wid41], Theorem VI. 19a) 


tp(m,X,y;C) = / due mu} ip(u, A, y; (). 


(5.7) 


Substituting (5.7) into the right hand side of (5.6) yields 


-^ 1 , 2 (Aq, A fc ) = -- y; 


1 w . „ T-(i+/t)m 2 -i m r°° 


1 — 1 crG{—1,1} 


m T \ 


due mu ijj(u,Xa,cr-,C)- 


In order to finish the verification of the identity for it 1,2 all that remains is to interchange 
the order of the summation in m and the integral in u, which, as required in view of (2.9), 
(2.10) and (2.12), leads to 

K$(X a ,X b ) = - £ / duF 3 (e~^(u,X a ,a-,0 

^ ^ r 1 11 J — OO 


with 


cr€{-l,l} 


“ T (\+^W-\m 

f 3 (z) = -j2 ———<- zr■ 

m— 1 


m T \ 


The application of Fubini’s theorem will be justified once we check that 

r-Ok+P )m 2 -\m 


/ OO 

du 

-°° m=1 


T y 4 


< 7 - 


m r 


-e mw 'ip(u,X,a;() 


< 00 . 


(5.8) 


We first observe that ip(s, A, cr; () = t~^ s2 ( s G(t s , A, cr), where G is a continuous function, 
whose arguments vary over compact domains (recall that s € iM). Therefore there is some 
C > 0 such that 

\ip(s, A, cr; C)| < Ct~P s2 \( s \. 

Let us write a = (i + /3) |log(r)| and b = (3 |log(r)|. Then we have, using the last inequality, 


, e bs 2 +soj _ 


\ip(u, A, cr; £)| < C f ds< 

J iR 

and thus the left hand side of (5.8) can be bounded by a constant multiple of 


^ e ~ uj2 / 4b i 


(5.9) 


/ CXD 00 

du 

-°° m= 1 


-am 2 — (1/4+cj) m—oj 2 /4b q 


, 1 


\fab J- 


f 


due( UJ+1 / A)2/4a - u2/4b < 00 


for some C' > 0, where the last integral is finite since 0 < b < a. 
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The argument for K\\ is similar. Using (5.5) and the definition of u ap in (4.7) we may 
write 

°° , ,, 2 


1 p 1 l Z 

K\i(X a ,X b ) = E - / %(-ir(l-r) 2 m r^+ 2 ^ 2 — -pL 

^i mJ c 0 ,i y 2 ~ 1 


< T l+m 2 \ t T l+mj 2 \ 

x V » oc ;°> (m ,A a ,y;CMm,A h 

( T y 2 'i T )oo \ T /y 2 'i T )oo y 

oc 

+ 9 E E (-!)-+-' (-aaO^+'sgn^W-am) 


m,m'=l o-,cr'g{ —1,1} 


r (i+/3)(m 2 +(m') 2 )-|(™+m') 


m r ! m(_! 


^(m, A a , a; Q^{m', X b , a'; C). 


Plugging in the Laplace transform formula (5.7) for i/>, we would like to interchange the 
integration in uj and oj' with the summation in m and m! as in the K\ 2 case. This leads to 

Ki,i(X a , Xb) = [ dujdu' — -f dyip(uj, X a , y; ()$((*}', X h ,hC)Fi(e~ u ~ u ' ,y) 

J r 2 m Jc 0 ,i y 

+ ^ E [ dw(L)'i>(u, X a ,(T;C)i>i^ > , X b ,a'-X)F 2 {e~ u ,e~ w ']a,a'), 


(T 1 <T , ^l{ — 1,1} ' 


with 


F ,t z „) - f' (- zra _ (P+’V;^ (P+’VlAr), 

“,SS s ' 2 - 1 WAHlVL 


T 2 (zi,z 2 ;ai,a 2 ) = E (~u) mi (~Z 2 ) m2 (-aia 2 ) miAm2+1 sgn(a 2 m 2 - aimi) 

mi,m 2 = l 

T ( | +P) (ml +m §)- \ (mi +m 2 ) 

X _ 

(mi) T ! (m 2 ) r ! 

which is exactly what we need. Hence all that remains is to justify the application of 
Fubini’s Theorem. To this end, in the case of the first summand it is enough to check that 



duj\ dui2 


1 /* 1 | y 2 — v 

ttw r dy „ i/>(o;i,A a ,y;C)^(u; 2 , Afc, ^;C) 

27ri 7c 0 i 3/ — 1 y 


E 

772 — 1 


3 — 1 + 072 ) 777 . 


r (i+2/3)m 2 -m 


(m r !) 


!}2 


< oo. 


(5.10) 


The y contour passes through singularities of the factor at y = ±1 and thus, as a 
principal value integral, it equals half the residues of the integrand at those points plus 
the same integral with the contour replaced by a circle Cut of radius slightly smaller than 
1. Now the residues of at y = ±1 equal ±1 so the residues of the integrand at these 

points equal ±t/>(u;i, A tt , ±1; £)^(u; 2 , Aft, ±1; £). On the other hand, the estimate (5.9) still 
holds (with a different constant) with a on the left hand side replaced by y, uniformly in 
y £ Coq U { — 1,1}. Therefore 


1 

27ri 



i + y 2 
y 2 -1 


V>(wi, A a ,y;C)^(w 2 , Aft, 



< Ce- uj2/Ab +Ce~ ul2/4b [ dy 
JC 0A 


i + y 2 
y 2 -1 


where, we recall, b = j5 |log(r)| (and later a = (-) + /3) |log(r)|). The last integral is finite, 
and thus this shows that the left hand side of (5.10) is bounded above by a constant multiple 
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of 



du i du 2 e -(^+“2 2 )/4fe 

m= 1 


e ~(uJl+uJ2+l)m 


g— 2am 2 

( m r !) 2 


< [ duJidu 2 e^ 1+U2+1)2/8a ~^ +UJ ^ /Ab , 

Jr 2 


which as before is finite thanks to the fact that a > b. This yields (5.10) and finishes the 
justification of the desired identity for the first summand of K\ \. The justification for the 
second summand is entirely analogous to the argument for Kf^ above, and thus we omit it. 
This concludes out proof of the equality between the kernels given in (2.12) and (5.3). □ 

6. Moment formulas for flat KPZ/stochastic heat equation 

In this section we will describe moment formulas for the solution of the KPZ equation 
with flat initial condition which are obtained in an analogous way as those for ASEP. 
Versions of these formulas already appeared in the physics literature in [LDC12], 

The one-dimensional Kardar-Parisi-Zhang (KPZ) “equation” is given by 

d t h = \d 2 h - \[(d x h) 2 — oo] + £. 

where £ is a space-time white noise. This SPDE is ill-posed as written but, at least on the 
torus, it can be made sense of by a renormalization procedure introduced by M. Hairer 
in [Hail3; Hail4], His solutions coincide with the Cole-Hopf solution obtained by setting 
h(t,x ) = — log Z(t,x), where Z is the unique solution to the (well-posed) stochastic heat 
equation (SHE) 

d t Z = \d 2 x Z + iZ. (6.1) 

The formulas which we will obtain will actually be for the moments E flat [Z(t, 0) m ] of the 
SHE with flat initial data, which means 

Z{ 0, x) = l or h(t,x) = 0. 

Of course, this means that we are getting formulas for the exponential moments E flat 
of flat KPZ. One can also think of this in terms of the solution of the delta Bose gas with 
flat initial data, which is the solution v(t]x) of the following system of equations, where 
we write H4 = {x €E : x\ < x 2 < • • • < Xk} (see [BC14; OQR14] for more details): 

(1) For x € W fc , 

dtv(t\x) = v(t\x), 

where the Laplacian acts on x. 

(2) For x on the boundary of Wk, with x a = x a +i, 

(dx a - d Xa+1 - l)v(t, x) = 0. 

(3) For x € Wk, 

lim v(t]x) = 1. 
t—>0 

It is widely accepted in the physics literature that if Z(t,x ) is a solution of the SHE, then 
v(t;x) = K[Z(t,xi) ■ ■ ■ Z(t,Xk)\ is a solution of the delta Bose gas. This fact is proved in 
[MFQR15], where it is also shown that there is at most one solution. 

In [OQR14] we obtained an explicit formula for the moments of the SHE started with 
half-flat initial condition by taking a weakly asymmetric limit of the half-flat ASEP moment 
formula (Theorem 2.4 above). The same formula was obtained in [OQR14] directly for the 
delta Bose gas with a slightly more general initial condition using the same method which 
led in that paper to the moment formulas for half-flat ASEP. 
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Proposition 6.1. For 6 > 0 define (with x appearing k times in the argument) 

k 

• 1 -—- i / 

dw I$(w ; n) 


vg- a (t-,x,...,x) = 2 


l\ ^ (2 'k\Y (a+mk 

1=0 ni,...,n/>l v 1 •R a + 1JK 1 

n\-\ -f -ri£=k 


( 6 . 2 ) 


with a > 0 and with Ig given by 


Ig(w,n) = 2 fc TT r(2w a ) ct[±n3-±n a +n a (wa-e) 2 ]+xn a (w a -d) 

-*■/ n a T(2w a + n a ) 

TT r(w a + w b + \{n a - n b ))T{w a + w b - \(n a - n b )) 

^ T(w a + w b - \{n a + n b ))T{w a + w b + \{n a + n b )) 

x K ~w b + \(n a - n b ))(w a -w b - \{n a - n b )) ^ 
(w a ~ w b ){w a -w b + \(n a + n b )) 

Then the solution v(t;x) of the delta Bose gas with the initial data in (3) replaced by the 
tilted half-flat initial condition lirri^o v(t; x) = [7,, ^ dXa l Xa >o coincides with x,... ,x) 

when x = (x,... ,x) (with i£l repeated k times). Moreover, if Z(t,x) is the solution of 
the SHE (6.1) with half-flat initial condition Z(0,x) = l x >o, then (denoting by the 
expectation with this initial condition) 

E h ' a [Z(t, x) k } = v$- Q (t; x,...,x). (6.4) 


This corresponds to Proposition 5.3 of [OQR14]. As in the ASEP case, we have used 
the Cauchy determinant identity (2.24) to expand the determinant which appears in that 
formula. 

Just as (6.4) was obtained in [OQR14] as a weakly asymmetric limit of (2.21), one can 
take a weakly asymmetric limit of (4.27) to derive a moment formula for the SHE with flat 
initial condition Z(0,x) = 1. Alternatively, one can take the x —>■ oo limit directly at the 
level of (6.2)-(6.4) to obtain such a formula. In fact, by statistical translation invariance 
of the white noise, Z{t,x) with half-flat initial condition Z(0,y) = l y >o has the same 
distribution as Z(t,0) with initial condition Z(0,y) = l y >- x , and thus computing the 
x —>• oo limit of (6.2) leads to the moments of Z(t, 0) with flat initial condition. This is the 
approach that we will take here (although we will not provide all the details). As for the 
case of ASEP, the calculations involved in the computation of this limit are quite involved, 
but they are completely analogous to those of Sections 3 and 4, as we describe next. 

Note that the dependence on x of the right hand side of (6.3) is only through factors of 
the form e XWarLa , which suggests that we should deform the w a contours to a region where 
lR(w a ) < 0. More precisely, we will shift the w a contours to — 5 + iR for some small 6 > 0. 
The contour deformation involves crossing many poles, and a careful analysis shows that 
the pole structure is entirely analogous to the one described in Section 3.1 for ASEP. Now 
the unpaired poles occur at w a = 0, coming from the factor T(2w a ), while the paired poles 
occur at w a = — w b in the case n a = n b , coming from Gamma factors in the numerator 
of the second line of (6.3). Crucially, as in the ASEP computation (though in a slightly 
simpler way), the dependence on x is gone in all the residues that one has to compute as 
the contours are deformed. Since the free variables will now be integrated on a contour 
—<5 + iR with 5 > 0, these factors will vanish as x —>• oo, and as a result the limiting formula 
will be written only in terms of unpaired and paired variables, analogously to what was 
done in Section 3 for ASEP. Instead of going through the whole argument again, we will 
only state the result. Introduce the index set 

A = ((n“, fl») € lix 

I a=l a= 1 
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Then 


E flat [ZM) m ] = 2 (6-5) 

fc =0 


with 


where 


-flat 

v k,m 


(*) = E 


1 




k u \2 k pk p ! 


H-2/cp—/c 


3 m*(^«Ap) — £ ttAt l n 

i 27T1 pp //'iTO'ifep J -- L 


P\2 


(nP,nP)€A- fcp 


a<fc p , b<k x 


(ng - - 4(zg) 

(ng + ng) 2 - 4:(za) 2 


K P ria J 


n 


(ng - <) 2 ~ 4(zq - ^ p ) 2 (ng - ng) 2 - 4(zg + zg) 2 
(ng + ng) 2 - 4(zg - z $) 2 (ng + ng) 2 - 4(zg + zg) 2 


l<a<b<kp 


— TT — 

}Alu -LX to U | 


x _ || _ e ^i(«) 3 -<) 

2 fc u 11 nil! 
a=l “ 


pj _L e i[iK) 3 -ing+ 2 nP(zP) 2 ] (_-q 

a=l l<a<6<fc u 


1 w£Ang K ~ I 


K + n b 

( 6 . 6 ) 


This formula should be compared with the moment formula for flat ASEP provided in 
Theorem 3.6 (in fact, (6.5) can be obtained as a weakly asymmetric limit of that formula). 
Note that in the above discussion the Za contours where deformed to —5 + iM, but here 
they have been replaced by iM; we may do this here (i.e. take <5 —>• 0) because the resulting 
integrand has no poles on z a G iM (this analogous to what was done for ASEP, see the 
paragraph following (3.20)). 

It is not hard to check that (6.6) corresponds to (108) in [LDC12] (with k u = M, 
k p = N and setting s = 0 in their formula, which comes from a scaling factor which we 
have not included at this point; note also that in [LDC12] the replacement t = 4A 3 has 
been performed). Their heuristic derivation essentially consists in studying the poles which 
arise from taking 6 —>• 0 at the same time as x —>• oo directly in (6.2) with a = 8. Although 
their procedure is different and not fully justified, the algebraic structure is similar. In 
terms of (6.2), they start by shifting all the w a variables by 9, which results in the contours 
turning into iM. As they take 0 —» 0, some of the factors in the integrand now pass through 
singularities, which give rise to residues which are in correspondence with the ones arising 
from our paired and unpaired poles. Then they argue that taking x —> oo at the same time 
as 6 —)• 0 results in only these types of terms remain in the limit, leading to the formula. 

The final step is to find a Pfaffian representation for the moment formula (6.5). This 
is done again in an entirely analogous fashion to what was done for ASEP in Sections 4.1 
and 4.2 (alternatively, and as we mentioned, one can take a weakly asymmetric limit of 
the moment formula in Proposition 2.3, which leads to the same expression). Let us skip 
the details and just write the result: 

171 ( 1 \k 00 r 

E^[Z(t,0r} = m\J2 [ -rr- E / dA Pf [it(A a , A fe ; m a , m 6 )] * 6=1 

k=0 * mi,...,mjfe=l ( R >o) fc 

raiH-fmpm 


(6.7) 
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with K{ Ai, A 2 ; mi, m 2 ) = 


Aa,l(Ai,A2;mi,m 2 ) 7(5,2 (Ai ,A 2 ;mi ,m 2 ) 
—A'l,2(A2,Ai;m2,mi) 7?2,2(Ai,A2;mi,m2) 


^i,i(Ai, A 2 ;m 1 ,m 2 ) = l mi=m2 (-l) mi+1 ^-r / 


and 

r(-2y) 


8vri7 iK 


+ 


1 


1 


1 


32 (mi — 1)! (m2 — 1)! 


y T(m + 2y) T(m — 2 y) 

X e [M m l + | m l27 2 ] t -l( m l- 2 27)Al-2(mi+2y)A 2 
(_l)”»iAm 2 S g n ( mi _ m2 ) glF2 ( m i+ m 2)* — j( mi Ai+ m 2 A2) 


^i,2(Ai,A2;mi,m2) = -i^^ T y I eTfe m ?*-i-^ 1 
^2,2(Ai,A 2 ;mi,m 2 ) = 5 sgn(A 2 - Ai). 

This formula should be compared with the the m-th term of the (divergent) series written 
for X) m >o ^ &at [Z(t,0) m ] (for a specific choice of C) in (147) of [LDC12] 5 . 

Appendix A. Formal GOE asymptotics for flat ASEP 


In this section we will provide a non-rigorous argument that shows that, under the 
scaling specified in (2.13), the right hand side of (2.11) recovers the GOE Tracy-Widom 
distribution. 

More precisely, we perform a formal asymptotic analysis of Pf[j — w ith 

C = —7 -~ t / 4 + tl/3r / 2 as t —> 00 and obtain Pf[J — K r \ L 2q 0o o))> with K r defined as in (2.15). 
This is the content of Section A.2, and should be interpreted as evidence for the conjecture 


lim P flat 

t—> OO 


( h(t/(q-p),t 2 / 3 x) 

\ t 1 / 3 



Pf[J — AT r ] L 2 {[0 oo)) . 


The issues discussed in Remark 2.2.vii preclude directly making this part of the asymptotics 
rigorous. However, even the formal critical point analysis has a complicated algebraic 
structure which needs to be explained. We will then show, rigorously, in Section A.3, that 
the right hand side recovers the Tracy-Widom GOE distribution: 


Pf[J - K r \ L 2 ([0)Oo)) = Egoe(t). 


A.l. Mellin-Barnes representation. There are two alternative ways to derive (formally) 
the t —y 00 asymptotics of our formula. The first one consists in starting directly with the 
formula given in Theorem 2.1, which is written in a form that lends itself readily to critical 
point analysis (this is thanks to the rewriting of (5.3) through (5.7), which is analogous to 
the physicists’ Airy trick, see the discussion in page 32 preceding the proof of the theorem). 
The difficulty with this approach is that it involves delicate asymptotics of the functions 
Fi, F 2 and F 3 (defined in (2.10)) which are as poorly behaved as exp T . 

A different argument relies on using a Mellin-Barnes representation (see [BC14], which 
introduced this idea in this setting) for the sums (in m and m') which appear in (5.3). 
The advantage of this approach is that the critical point analysis turns out to be relatively 
simple and one is not faced with the asymptotics of the functions F\, F 2 and F%. The 

■'’There are two main differences between (6.7) and the n-th term of the formula written in [LDC12], 
besides some simple scaling and renaming of variables. First, in their formula they have applied the Airy 
trick (5.4) in order to rewrite expressions involving e® m t in terms of integrals of Airy functions. This can 
be done for (6.7) without any trouble (we have not done it because it leads to more complicated formulas). 
Second, one can check that obtaining their formula from ours involves (formally) integrating in Ai and 
A 2 , integrating 1 in Ai and differentiating it in A 2 , and differentiating K. 2,2 in both Ai and A 2 . This can 
be justified (formally) using Lemma B.3 (twice). This difference stems from the delta Bose gas analog of 
the step performed at the ASEP level after (4.19), which is not done in [LDC12]. 
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difficulty, on the other hand, is that implementing the Mellin-Barnes representation turns 
out to be much harder than usual, due to the existence of a large number of additional 
poles (on top of the ones poles at the integers which are inherent in the representation) 
which have to be accounted for. The fact that all these additional poles end up cancelling 
is far from trivial and serves as another indication of the remarkable structure behind the 
flat ASEP formulas. Partly because of this reason, this is the approach we will follow 
presently. 

In order to derive the Mellin-Barnes representation of the kernel A" flat, £ in (5.3) we will 
need to use the following generalization of the Mellin-Barnes representation formula which 
appears e.g. in [BC14]: 


Lemma A.l. Let C\ 2 ... be a negatively oriented contour enclosing all positive integers 
(e.g. Ci t 2 ,... = \ + iM oriented with increasing imaginary part). Let g be a meromorphic 
function and let A be the set of all poles of g lying to the right of C\ 2 ...- Assume that 
A n Z>i = 0. Then for ( £ C \ M>o with |C| < 1 we have 


OO 1 p 

Yg(n)( n = - / 


ds 


IT 


sin(— ns) 


(-c)*sw - E 


7r 


ipGA 


sin(— irw) 


(-erRes g(r s 


provided that the left hand side converges and that there exist closed contours Ck , k € N 
enclosing the positive integers from 1 to k and such that the integral of the integrand on 
the right hand side over the symmetric difference of C 1 . 2 ,... and Ck goes to zero as k —>• 00 . 


The proof is elementary and only uses the fact that 7 r/ sin(— tts) has poles at each s = 
l:£Z with residue equal to (—l) fc+1 . The (trivial) novelty of this formula is the appearance 
of a second family of poles at A. 

In order to apply this result to our formula we need to find a different expression for 
which as written contains factors of the form 1 ma=mb , (—l) m “ Am6 and (—l) m “ 
(appearing inside the definition of u ap (y,m a )) which are not suitable for a representation 
in terms of contour integrals. The following result can be proved in a similar way to Lemma 
3.1: 


Lemma A.2. Given any m\ m 2 € Z with the same sign, 

lim r AAl - + fLip ~ mi) + "\ - mi + 2 " = (—l ) miAm2 sgn(m2—nii) 

T(i(mi + m 2 ) + 7/)r(— i(mi + m 2 ) + rf) m\ + m 2 + 2 rj 


(One can choose slightly different versions of the left hand side, but this form will turn 
out to be convenient later). As a consequence, if we define 


h f? (o‘i, <t 2 ; mi,m 2 ) 


r(i(mi - m 2 )+ 7 ?)T(i(m 2 - mi)+ 77 ) m 2 - mi + 2 r/ 

(79 ——j—---II (7f — (J 2 

< T(i(mi + m 2 ) + 7 /)r(—±(mi + m 2 ) + rf) mi + m 2 + 2 77 
o 2 if cr!^cr 2 , 


then we may replace every factor of the form (— cr a ab) maArnb sgn (crbmb — cr a m a ) in our 
formula by lim^o ^(o'a, m a , m&). Something similar can be done for the product of 
1 m a =m b and the factor (-l) m “ appearing in u ap (77 a ,m a ), replacing it by lim^ 0 s r? (m a , m b ) 
with 


s,(mi,m 2 ) 


sin(^7r(2mi + 1)) 


sin(7r(mi — m 2 + rf)) 
7 r(mi - m 2 + 77 ) 


As we will see in the next result, the 77 —>• 0 limit can be taken outside the sums and 
integrals. We state it at the level of the formula for the exp T (x;C) transform: 








A PFAFFIAN REPRESENTATION FOR FLAT ASEP 


40 


Proposition A.3. For every £,£ € C with |£| < t 1 / 4 or |£| = r 1 / 4 and |£| < t 1 / 4 we have 


E : 


flat 


exp T (CTl^o) ;e )j = ^ E( _ 1)fc _ £ 


rj—>0 ■ 

k =0 


d\ 


m i ,...,771/2 


= i •7 (R>o)' £ 


X ^E a ^a T -i(E a m a )^(E a m Q ) 2 -E a m a p f A^’^A^A^m^mb) (A.l) 

J a, 6=1 

where A flat,? ? is tde 2x2 skew-symmetric matrix kernel defined by 

K^f\X a ,X b -m a ,m b )= 5v (m a ,m b )--f dy r^(X a , y, m a )t>(A 6 , ±, m a )u ap (y, m a ) 

vnJco.i 

+ 77 ^ f)r?(o-, cr'; m a , m b )T^ m * +m ^ d(A a , a, m a )d(A fe , a', m b )u u (cr, m a )u u (o-', m b ), 


cr,ir'e{—1,1} 


r.- ftfll . !j 7^flat 

^ 1,2 — -^ 1,2 ) 


and 


flat,77 _ t -flat 


A 22 ’' = 


2,2 1 


where 


u ap (y,m) = (—l) m u ap (y, m) = r 


(A.2) 


and u u and 0 were defined respectively in (4.7) and (4.23). 


Proof. Consider the kernel A flat, ° defined as lim^o A flat,T h Then it is straightforward to 
check from Theorem 5.1, using Lemma A.2 and the fact that s v (mi, m 2 ) —> (—l) mi 1 m 1 =m 2 
as 7 ? —)■ 0 for mi, m 2 € Z, that 


E 


flat 


exp r (CTf'-«'»>;£)l £ 


k =0 


,.^ fe =l ( 27ri) • / ( R >o) fc 


dA 


X ^E a m ar -i(E a m a )^(E a m a ) 2 -E a ma p f [A flat ’°(A a , A fe ; m a , m b ) 1 * 


- a,6=1 


To justify taking the 77 —>• 0 limit outside use the Dominated Convergence Theorem together 
with an argument similar to the one used in the proof of Theorem 5.1 to show that the 
series on the right hand side of (A.l) is absolutely summable. □ 


Note that, since s v (m a ,m b ) is replacing the factor (— l) ma l ma=mb , the first term of 
K^'f ,r> in Proposition A.3 can be rewritten as 

1 [ f °° 1 2 

s r] (m a ,m b )C ma ~ mb — / dy dv ( 2 ma T 2 m "-v(Xa, y, m a )t>(X b , l/y, m a )u ap (y, m a ) 

714 Jc 0 ,i do 

(as can be checked directly in the proof). Using this replacement, setting £ = r 1 / 4 in 
Proposition A.3 and bringing the sums in m a inside the Pfafhan as in Section 5 yields the 
following formula: 


E : 


flat 


exp. 


_(£ T |fe(t.O)) 


= lim Pf 

J-K 71 

77 —>-0 

■ 


(A.3) 
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with (we omit here the dependence of the kernel on £) 


1 /’ 

K i,i( x , x ') = Y s v {m,m')C m ~ m '— f dy y, m)t>(A', ± m)u ap {y, m) 




m+m 1 j (m 2 +(m') 2 ) 


m,m’=1 o-,o-'e{—1,1} 


x d(A, cr, m)d(A', cr', m')u u (cr, m)u u (cr', m'), 


^i, 2 ( a , A') = -i Y aC m r4 m2 d(A,(7,m)u u ((7,m), 

m=l crG{—1,1} 

^ 2 (A,A') = isgn(A / -A). 


Up to here we have proceeded rigorously. In what follows we will proceed with a formal 
asymptotic analysis. The first step, which we undertake in the rest of this subsection, is to 
use a Mellin-Barnes representation (Lemma A.l, however, without verifying the necessary 
decay conditions) to argue formally that the right hand side of (A.3) yields 


Pf 


J-K° 


(A.4) 


with 

Ati(A,A') 


/ii° 2 (A,A') 

KlAW) 


1 


f . 7T 2 

/ dsds T-. -——- tS 0 (s,s) 

j(c+ iR )2 sin(- 7 rs) sm(— 7 rs') 


( 27ri ) 2 y ( c +*) 2 

x C s ~ s '-^i dy y, s)d(A', ± s)u ap (y, s) 

ni Jc 0A y 

1 f 7 T^ 

+ -x 2 / dsds 1 — -——-- V fio(o-,cr';s,s') 

2(27Tl)- 7( c+iR )2 sm(-7rs) sm(-7TS / ) aa ,^r J _ l 


■— f 

4vn .7 c+i 


x (-cr +s r 


s+s' ±(s 2 + (s') 2 ) 


ds 

c+iR sm (“ 


6{-l,l} 

o(A, £ 7 , s)d( A', cr', s')u u (cr, s)u u (cr', s'), 

14 


7 - 7 Y a (~0 S T 4S 0(A,CJ,s)u u (£J,s), 


5 sgn(A' - A), 


(A.5) 


where So(s, s') = lirri^os r? (s, s') and flo(cr, cr'; s, s') = lim^o f)r/(cr, cr'; s, s'). Note that this 
formula involves using first the Mellin-Barnes representation and then computing the limit 
r] -A 0. As we will see below, it is only in the limit rj —7 0 that the additional poles which 
arise in the Mellin-Barnes representation cancel. 

In view of (A.3), and taking the 77 —>- 0 limit back inside the Pfaffian, our goal is to argue 
that lim^o is given by K °. Consider first applying the Mellin-Barnes representation 
to the kernel K y 2 . One can check that in this case the only poles of the integrand (in s) 
are the ones occurring at s € Z>i coming from 7 r/ sin(—7rs), and thus the Mellin-Barnes 
representation can be applied without additional difficulties (in other words the set A in 
Lemma A.l is empty in this case). Doing this and taking 77 —>• 0 yields K® 2 - Note also that 
K 2 2 does not depend on rj, and in fact it equals K 2 2 , so the equality is direct in this case. 

It remains to handle which is where the additional poles in the Mellin-Barnes 

representation arise. Write for convenience 


*?,i =L\ + L\ 
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where L/ and L r / correspond to each of the two terms appearing in the definition of K r / j. 
For L r { one checks again that the only poles occur at s, s' € Z>i so the Mellin-Barnes 
representation can be applied as before, yielding 


L i = 7 ^—Y 2 dsds' — --- j 7 S v (s,s') 

( 27ri ) J(c+ iM )2 sm(-7rs) SUl(—7TS'J 


_n—-'ll 


X (-0 


-J dyC T 2 S o(A,y,s)t)(A , s)u ap (y,s) (A.6) 

m JCo,i 


for c € (0,1). Taking r) —> 0 will yield the desired representation. 


Now write L 2 as 


L 2 = 5 X X] a') 

m,m'= 1 ct,(t'G{—1,1} 


(A.7) 


with = r A m +( m ') )t)(A, cr, m)t>(A', cr', m')u u (cr, m)u u (cr', m'). Our goal is to 

get 


lim 

97^0 z 


2(2vri) 


ds ds' 


ir 2 (-() s+ s' 


'(c+iK) 2 


sin(—7rs) sin(— ns') 


X i)o{s,s';a,a')I{s,s'-,a,a') 


( 7 ,( 7 ' — 1,1} 


This case is more involved, due to the factor ^(s, s' ; cr, a'). One can check that I(s, s'; cr, cr') 
is analytic away from s, s' € Z>i, and in fact the same holds for (^(s, s'; cr, cr') when a / cr'. 
Now consider the case a = a'. Let us first fix m and consider the effect of applying Lemma 
A.l to the m' sum in (A.7). We need to analyze the poles (in s') of 

T(i(m-s') + ? 7 )r(^(s'-m) + ??) s'- m + 2 r] 

n{m,s ,a,a) a r (i + + ^ r (_ 1 ( m + a /) + v ) m + a ' + 2rj ' 

The numerator has singularities when s' = m+2£+2r/ and s' = m—2£—2r] for £ € Z>o- The 
first type of singularity is removable, because at that point the factor 1 /T(— \{m + s') + rj) 
evaluates to 1/T (—m — £) = 0. The second type of singularity is, on the other hand, a pole, 

with residue 2 o ^~I r(rnS)r(e-m+ 2 r)) ’ an< ^ w ^i c ^ 1 belongs to the deformation region 

(namely {z € C : 9 i(z) > -(}) only when £ < m/2 (assuming that y is small). In view of 
this, Lemma A.l suggests that 


L\ = 


X / ds ' — A X C m (-C) s I» 7 ( cr i cr, ; m >' s, )7(o-,o-';m,s') 

^£1,1} 

(~ll T(^ + 2 rj) 1-2 77 ^ ^ m _ 2e _ 2ri 


00 f m /2j 

EE E ■ 

m =1 £=0 erg{—1,1} 


H T(m — ^)T(^ — m + 2?y) ^ — m 


7T 


sin(7r(2ry + 21! — m)) 


/(cr, cr'; m, m — 2 t — 2 77 ) 


— 7^ - 

iJ 2ll -^ 2 , 2 - 


Here we are choosing the contour C\ 2 in Lemma A.l to be c + iM. Now we need to apply 
the Mellin-Barnes representation to the remaining sum (in m) appearing in L 1 / 1; regarding 
s' as fixed. As before most of the factors making up the integrand are analytic, except for 
t) v (s,s';a,a') when cr = a', in which case it reads 

T(±(s - s')+ 7?)r(±(s'- s)+ry) s '-s + 2ri 
S,a,a) CT r(i(s + s') + r?)T(-i(s + s') + 77 ) s + s' + 2 y 

The numerator has singularities when s = s' — 2 £ — 2 y and s = s' + 2 £ + 2 y for l € Z>o- 
Note that the singularities of the first type are never in the deformation region, so we are 
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only left with the second type of singularities, which lie in the deformation region for all 
£ > 0, with residue 2a^- r ( s / + £+ 2 ^)F(-^-£) s'+e+2 V - Reasoning as before we get 

r-2( /-\S + S' 


, , 1 f Tr 2 ( _t\ S+S 

L » = ww dsis ' *■(-») It -,,-) 1 a ' w> ’ s,; * <7 ' ) 
"5; J-r.,2^/+u 


D' 


r(< + 2 t?) 


£>0 ctG {— 1 , 1 } 


t\ t( s ’+ £ + 2r])T{-s'-£)s' + £ + 2rj 


2(_Q2s'+2e+2r) 


- 7 - T~i -7 , -tt- 7 -— I (s' + 21 + 277, s'; < 7 , a) 

sm(— 7 r(s' + 2 £ + 2 r/)) sm(— 7 rs') 


:= L\ z - V 


2,4- 


In terms of the above kernels we have K\ 2 = L r ( + L 2 3 — L 2 2 — T 2 4 - Observe now, using 
(A. 6 ) and the last equation, that the limit as 77 —>• 0 of L 4 + L 2 3 yields exactly K ®, (given 
in (A.5)), so all that remains to show is that L 2 2 + L 2 4 —> 0 as 77 —> 0. For the first term, 
noting that the zero of the sine in the denominator cancels with a zero coming from one 
of the Gamma functions, so that si a (n( 2 ' n + 2 £-m))r(e-m.+ 2 r 1 ) —* ~ we have: 


00 L m / 2 J 

;™ l 2,2 = -E E E ^-irc 

m= 1 £=0 <re{-l,l} 


m/-2m—2l 


I (a, < 7 ; rn, m — 2 £) 


(A. 8 ) 


EE E 0^_-j^m^2m+2£ J(^j cr" 777 —F 2 £ 

£=0 m =0 cre{—1,1} 


m) 


where, for convenience, in the second equality we have added the term with £ = m = 0 
which is 0 anyway because u u (< 7 ,0) (which appears in I(a, cr; 0,0)) is so. On the other hand 
we have 


lim l 2 4 = y - [ i 

2 ’ 4 &-r.fi a 


ds / (7 


-it 


r(«) 


£>0 crG{—1,1} 


£! r(s' + ^)r(-s'-£)s' + £ 
7 r 2 ( _ C) 2 ,'+ 2 £ 


X 


sin — 


7 -—--rr- 7 ---/(s / + 2 £, s 7 ; cr, a). 

(— 7 r(s' + 2 £)) sm(— ns') 


The s' integral can be computed in terms of the residues of the integrand for ^. 

There is a double zero in the denominator coming from the sine factors when s' = m! € Z>i, 
but one of them is canceled by the zero of 1 /T(— s' — £) at these points, resulting in a simple 
pole, with Res s ' =m 7 r 2 /[sin(— n(s' + 2£)) sin(— 7rs')r(— s' — £)\ = + £)\. We get 

lim L 2)4 = E E E 7 (-l) m '( 2 m ' + 2 £ rJ« m ' + 2 f ’ 2 R m ') 2 )/(a, a; m' + 2 f, m'). 

V «>0m'>l(rS{-l,l} 

As before we may add the term with £ = nn! = 0 , and now comparing with (A. 8 ) we see 
that lim^o(E 2 2 + L 2 4 ) = 0 as desired, which finishes our derivation of (A.4). 

A.2. Computation of the limit. We are finally in position to compute the t —» oo 
asymptotics of the flat ASEP distribution function. We take £ = —t~ 4 t+ 2 i /3r and, in 
view (A.4), we need to compute the limit of Pf[J — K°], with K° given in (A.5). We will 
only provide a formal critical point derivation of the limit. To that end we study first 
the factors in K°(s,s') which depend on t. They come from the products of the form 
r ( ~3*+^ 1 / 3 r ) s t)(A, y, s), and are given by 

l 


exp 1 1 


1 +t -s/2 


y 


1 + rs/2 y - 3 log(r)s + \t 1/3 rs log(r)) 
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(this factor appears twice in A' 3 l (A. A'; s, s'), once with respect to s and once with respect to 
s ', and once in K® 2 { A, A'; s, s')). Consider now the function f(s, y) = 1 +T ~ s / 2 y ~ 1 + ^s/ 2 y ~ 
i log(r)s. One checks that / has a critical point at (0,1) and, moreover, that the Hessian 
of this function vanishes at this point. This suggests that we should use a f -1 / 3 scaling 
around these points. Explicitly, we will rescale as follows: 

y = 1 + t~ 1/3 y , 


<. _ -1 /—1/3— 


and 


A = t 1 / 3 (A-2r) 


1 s/2 

(and similarly for X' and s'). The A scaling is so that A |/^ a / 2 ^ is order 1 (the shift by 2 r is 

for convenience). The A change of variables produces a t 1 ' 3 in front of the kernel K°, and 
using (B.4) this can be rewritten as 


t 1 / 3 iAT°(A, A') = 
With this scaling we have 


t 2 / 3 K° A (\,X) W 3 *r° 2 (A,A') 
A) K°J\ A') 


(A.9) 


t 

i-W 2 


1 + T~ S / 2 


hj 27 , - i+X 3 y ~ 2 lo s( T ) s l + \ rtl/i l °s(r)s « ^s 3 + ^ 


2/ 


16^ 


i s, 2 V ~2 (s- 2 y)t~ 1/3 , T ~ a / 2 lpi-dy w ±dy, t~V 3 . A ' daw-Urfs, 

l-j— j~ s /^y 4 v. > j/^ — 1 a y sm(—7rs) s ’ 

sin(7r(si—S 2 )) 

1-1 ^ 5 v — 2 v 


ein/'ZLfOci 4- 1 'A gm^qsi-sgjj ^ i r(§(si-s 2 ))r(lQ 2 -si)) S2 _ S1 ^ si+s 2 

V 2 1 “' 1 ' )) 7 r(si-s 2 ) ~ ’ r(|( Sl + S 2 ))r(-i(5i+ S2 )) S1+S2 ~ S2-S1 


and, furthermore, using Lemma A.l in [OQR14], 


(-r s/2 y,r) c 


and 


(O+Vu), 


1. 


( -T 3 / 2 y;r) oo ~ ^ (ry' 2 ;T) a 

Using these asymptotics and the definition of n in (4.23), we get 

li(A, y, s) » it _1/3 (s - 2 g 0 eife s3 +^-i( s '- 2 y)A+|yr 
= -t- 1 / 3 a A ei52 S3+ ^ s ^ 2 -^(*- 2 ^ + 5^. 

When y is replaced by 1 /y the above asymptotics hold with y replaced by —y. At the same 
time, if y = 1 then the asymptotics for 0 hold with y = 0 while if y = —1 the right hand 
side should be replaced by 0. The above asymptotics also give, in view of the definition of 
u u in (4.7) and of u ap in (A.2), 

u u (±l, s) « 1 and u ap (y, s)dy ~ |dy. 

Using all this in (A.5), choosing the s and s' contours to be ^ t _1 / 3 + iM, and keeping (A.9) 
in mind we deduce that 

(A.10) 


lim K° = K 

£—>■00 


with 


ATh(A,a') » ardr/i/ dyj^-f dSdS 1 Xle^+^-^+^+hi^y 
M f A A m lR y ^Ja +m) 2 ss y 


+ %A'iw/, 

J(h+a 


K^( A,A # )«-^ dSe ife 53 -^, 

/ 1 i ;td> 


2ss / (s'—s) 


AT 2 , 2 (a, A ) = i sgn(A - A). 

What remains is to rewrite AT in a more convenient way (this part bears some similarities 
with Appendix K of [LDC12]). Let us write Ad,i = A"^ + A'x.i- Note that K ] Li depends 
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on s' only through the integral 2 ^j;fi +iR ds' ^ which, as a principal value integral, equals 
-7. Then, using the definition of the Airy function, 

Ai(z) = J-f due^~ uz (A. 11 ) 

2m J c+m 

for c > 0 , we get (removing the tildes) 


i^Ai,A 2 ) = d Xl d X 2 X-f dy+'f ds 1 ^ e M+^y 2 -Usi-iy)M-l(s 1 + 2 y)x 2 
Jm. Jc+m 

= d A^i{ dy f°° du Wi f ds^eM+^y 2 -^-^- 1 ^ 

JiR JO J c+iR 

p p 00 

= d Xl d X2 ± [ l dy / ^iAi(2 2 / 3 (i(A 1 + A 2 ) + y 2 )+u)e i ( A2 - A ^ 

Jr Jo 

p poo 

= dx 2 ^hfdy / du i Ai / (2 2//3 (i(A 1 + A 2 ) + y 2 ) + u)e^ X2 ~ Xl ^ y 
Jr Jo 

p poo 

- d x 2 h T dy / du Ai(2 2/3 (i(Ai + A 2 ) + y 2 ) + u)d 
Jr Jo 

= -^a 2 ^- [ dy -f dy Ai(2 2 / 3 (i(A! + A 2 ) +y 2 ))e" 

J Ai-a 2 Jr 


(A 2 — \i)y 


p poo 

- 9 X 2 2 hi dy j du Ai(2 2 / 3 (i(Ai + A 2 ) + y 2 + 2 ~ 2 /\))e^- x 9 y _ 
Jr Jo 

As in Appendix J of [LDC12], we will use the following identity, proved in [VSI97]: 

/ OO 

dy Ai(2 2 / 3 (t / 2 + i(a + 6 )))^"^ = Ai(a) Ai(b). 

-OO 


It implies that 


K\ i(Ai, A 2 ) = ~d X2 \ [ dr/Ai(i(Ai + A 2 + ?y)) Ai(i(Ai + A 2 -r/)) 
J Ai-a 2 

roo 


poo 

- 9 X 2 ^J 3 J du Ai(Ai + 2 " 2 / 3 u) Ai(A 2 + 2 ~ 2 ^u) 

\ [ dyA.i{rj) Ai(Ai + A 2 - 77 ) - ~ Ai(Ai) Ai(A 2 ) - d X 2 K^{\ x , A 2 ), 

" J Ai 

(A.12) 

where K\\ was defined in (2.16). 

Now we turn to K 2 15 which is given by 


kK++) = a Xl d x ^l ds , 

J (c+iR) 2 

= e Xl e x ,^ f ds , / *2 SI m^eife(*S+4)-i<..A,+. 3 A 2) _ (Al , , a 2 ) 

v ' J2c+iR J c+iR ^ ; 


where we have shifted the s 2 contour to 2c + iR for convenience and where (Ai <—> A 2 ) 
denotes the same as the previous expression with Ai and A 2 interchanged. The first term 
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on the right hand side can be rewritten as 

[ ds i / ds 2 [°° du i [°° du 2 e TS3( s ?+*iM(»i*i+«A2)-(-i-«)ui-MU2 

J 2c+iR J c+iR J 0 JO 

= -<9a+a 2 |t+w / dsi / ds 2 /°°dPi r dv 2 e^ 2 ^ +s ^-^ SlXl+S2X ^- SlVl - S2V2 
J 2c+iR J c+iR J 0 J—v\ 


POO POO 

= —8d\ 1 d\ 2 / dv i / dv 2 Ai(Ai + 4vi)Ai(A 2 + 4v 2 ) 

JO J—v\ 

POO POO 

= \ / dv Ai'(Ai + v) Ai(A 2 — v) = | / dv Ai'(v) Ai(Ai + A 2 — v), 

J o " y Al 

where we have used again (A.11). Subtracting the term with Ai and A 2 flipped we get 

Afi(Ai,A 2 ) = \ f du Ai^cu)Ai(Ai + A 2 — uj). 

J \i 

Note that this term cancels the first term on the right hand side of (A. 12), and thus we 
obtain 


A'i,i(Ai,A 2 ) = — <9 A2 ATAi(Ai, A 2 ) - ^ Ai(Ai) Ai(A 2 ) = ^(<9 Al — 5 Aa )AAi(Ai, A 2 ), 
where the second equality follows by integration by parts. 

K\ 2 is simpler to obtain: using (A.11) one more time, 


Au, 2 (Ai, A 2 ) 


i 

167ri 


ds e M 2 


T i 7 s 3 -isAi 


' c+iR 


-i Ai(Ai). 


This, together with (A. 13), gives 
A(Ai,A 2 ) = Is 


Hd A, 


— (9 A2 )AAi(Ai, A 2 ) 
A Ai(A 2 ) 


-A Ai(Ai) 
A S gn(A 2 - Ai) 


(A.13) 


In view of (A.4), (A.10) and the definition of K r in (2.15) (shifting Ai i-A Ai + r, A 2 eA A 2 +r 
in the Pfaffian as well), this finishes our formal asymptotic analysis of Pf[j — / A flat >CJ, 
yielding Pf [J — AT r ], 


A.3. Fredholm PfafBan formula for GOE. All that remains is to complete our formal 
derivation of the GOE asymptotics for flat ASEP is to verify that the limiting Fredholm 
Pfaffian is indeed a formula for the Tracy-Widom GOE distribution. This fact can be 
proved rigorously: 


Proposition A.4. 


with 

A r (A l5 A 2 ) 


Agoe(t) = Pf[J - A+2 ([0iOo)) , 

A(<9ai “ <9 A2 )Ar A i(Ai + r, \ 2 + r) Ai(Ai + r) 
AAi(A 2 + r) A sgn (A 2 -Ai) 


(A. 14) 


A very similar formula for Agoe appears in [Fer04]. Verifying that the Fredholm Pfaffian 
on the right hand side of (A. 14) is convergent is not hard. In fact, it is enough to note 
that the kernel K r is uniformly bounded for Ai, A 2 > 0, and then expand the Pfaffian as a 
Fredholm series and use the identity Pf[A] 2 = det(A) together with Hadamard’s bound. 

In view of the discussion about Fredholm Pfaffians contained in Appendix B.3 (see in 
particular (B.17)), it would be nice if one could find a symplectic 2 x 2-matrix kernel M 
so that M J K r M is trace class. Unfortunately, it does not seem like such a kernel exists. 
More precisely, at least in the Fredholm determinant case one is usually led to consider 
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multiplication operators for conjugation (see (B.10)). In the symplectic case such an op¬ 
erator would be of the form M(^)( x) = f° r some (non-vanishing) function 

(/). But it is proved in page 82 of [DG09] that there is no choice of </> for which M T K r M 
is trace class (the underlying problem is the singularity at the diagonal of sgn(A 2 — Ai)). 
As mentioned at the end of Appendix B.3, we can still use the relation between Fredholm 
Pfaffians and Fredholm determinants since we know that the Fredholm Pfafhan in (A. 14) 
and the associated Fredholm determinant (see below) define absolutely convergent series. 
This is what we will do in the proof, which will then be devoted to showing that the re¬ 
sulting Fredholm determinant yields Fgoe- The proof follows the arguments of Section 7 
of [LDC12] relatively closely. 


Proof of Proposition A.f. For convenience will omit the subscript r in K r during this proof. 
We will also omit the subscripts from Fredholm determinants and Pfaffians, which are 
always computed on L 2 ([0,oo)) or L 2 ([0,oo)) <8> L 2 ([0, oo)). 

We have shown already that the Fredholm Pfafhan series for Pf[J — K] is absolutely 
convergent. An identical argument shows that det[J + JK ] satisfies the same, and thus 
Proposition B.4 implies that 

Pf[J — K] 2 = det [I + JK] = det 

In view of (A. 13), let us write 

ATpi = Kf tl + K\ a 

with 


I-KJ 2 K, 2,2 
-ATi I-K,2 


(A.15) 


Kf i (Ai, A 2 ) — —d\ 2 Kw(\\ + r, A 2 + r) and K 2 ^(\\. A 2 ) — —^ Ai(Ai + r) Ai(A 2 + r). 

Note that K\ 1 is a symmetric, rank-one kernel. We claim that K\\ can be replaced by 
K° j on the right hand side of (A.15): 


Pf [J - K] 2 


det 


I-Kj 2 A' 2,2 

-K{ A I—K 1} 2 


(A.16) 


It is enough to check this at the level of the finite dimensional determinants which appear 
in Fredholm determinant series (B.7) for det[7 + JK], in which case it follows from the 
following general fact: if A, B , C, U are n x n (real) matrices with B and C skew-symmetric 
and U symmetric and rank-one, then the matrices 


A B 
C + U A T 


and 


A B 
C A T 


(A.17) 


have the same eigenvalues (and eigenvectors) To see this, suppose that (iq, V 2 ) is an eigenvec¬ 


tor of 


A B 
C+U A T 


, where Vi € M n . We claim that Uv 1 = 0, which implies that (ui, U 2 ) is a ls° 

an eigenvector of ^ T ] with the same eigenvalue, say A. In fact, we have Av\ +Bv -2 = Xv\ 
and ( C+U)v\+A T V 2 = Xv 2 - Now test the first equation on the right with V 2 and the second 
on the left with v\ to get (Av\,V 2 )+(Bv 2 ,V 2 ) = A(vi, V 2 ) and (v\, {C+U)v\)+{v\, A t Tv 2 ) = 
A(ui,U 2 ). By skew-symmetry we have (Bv 2 ,^ 2 ) = (vi,Cvi) = 0, so we conclude from the 
two equations that (ui, Uv 1 ) = 0. Now a symmetric rank one matrix must be a multiple of 
a projection onto a vector w and we conclude that {v\,w} = 0 which further implies that 
Uv 1 = 0. This proves our claim that the matrices in (A.17) have the same eigenvalues, and 
thus (A.16). 

Since K\^ is rank-one, the kernel inside the determinant on the right hand side of (A.16) 

1 A'2,2' 


— K a 
^ 1,1 


Now recall that if B is a rank-one kernel then 


is a rank-two perturbation of 

det [I + A + B] = det [7 + A] (l + tr[(I + A)~ l B]) 


(A.18) 
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(this will be used later on), as long as I + A is invertible. A similar formula holds for 
rank-two perturbations, but in our case the skew-symmetry of Kf 1 and K 2>2 yields the 
nicer formula 


Pf [J - K ] 2 = det 


1 + 


0 

- K ? 


p2,2 

0 


Kj,2 


0 

0 Ki 2 


= det[l + (1-tr [{I + K^]) 2 . (A.19) 

Of course, a necessary condition for this identity to hold is that I + K'{ x K 22 be invertible, 
but this is true thanks to (A.22). The identity follows by approximation from the following 
matrix identity, which is (167) in [LDC12]: if A and P are skew-symmetric n x n matrices, 
R = uv T for some «,v£R", Q = [ ^ ] and P = [ ^ f ], then 

det(P + Q) = det(P)(l + (it, (/ - AB)~ l v)) 2 , (A.20) 

as long as I — AB is invertible. Since this matrix identity is not completely obvious, and 
no proof is given in [LDC12], let us pause for a moment to prove it. 

Note first that, since I — AB is invertible by assumption (and then, taking transpose and 
by skew-symmetry of A and B, so is I—BA), we have P -1 = (y _, [ J A ~P ]. 

On the other hand, we have det(P + Q) = det(P) det(7 + P~ 1 Q), so it suffices to show that 
det (I + P~ l Q) = (1 + (u,(I — AB)~ 1 v)) 2 . Using the formula for P~ 1 we get (P~ 1 Q)(wl) = 

(v,wi)(I—BA)~ 1 u 
(u,W2)(I—AB)~ 1 v 


((v’wi) fi-AB)~ 1 Au )' N° w since A and B are skew-symmetric, we have 

(taking transpose) that (u, (/ — AB)~ l Au) = —(u,A(I — BA)~ l u). But (I — AB)~ 1 A = 
A(I — PA) -1 , as can be checked easily by multiplying both sides by I — AB, so the 
previous identity implies that (u, (/ — AB)~ 1 Au) = 0. Now it is easy to see that P~~ 1 Q 


is rank-two, with eigenfunctions 


(I—BA)~ 1 u 

(I—AB)~ 1 Au 


and 


(. I-BA)~ 1 Bv 
(I—AB)~ 1 v 


, both with eigenvalue 


(u, (/ — AB) 1 v). This implies that det(I + P 1 Q) = (1 + (u, (I — AB) 1 v)) 2 as desired, 
and (A.20) follows. 

Going back to (A.19), observe that 

roo 

7Ci“iK 2 ,2(Ai,A 2 ) = -i / d£%RT A i(Ai+ r,£ + r)sgn(A 2 -0 

J 0 

= -K a i(Ai + r, A 2 + r) + \K a i(Ai + r, r). 

Writing 

Pr(Ai, A 2 ) = Ai(Ai + A 2 + r) and <5 0 = do ® 1 


(here 1 is the function which is identically equal to 1), observing that K Ai (- + r, - + r) = B 2 
and using the last formula we may write 

K? :1 K 2 ,2 = -B? + ±B% and K 1>2 = -\B r 5 0 . 

The kernel \B 2 8 0 is rank-one, and thus by (A.18) we have 

det[J + Kl x K 2)2 ] = det [I - P 2 ]( 1 + itr[(J - P 2 ) - 1 P 2 j 0 ]). (A. 21 ) 


On the other hand, using again the fact that Kf l K 22 is a rank-one perturbation of — B 2 , 
the Sherman-Morrison formula gives 


(I + k° xa k 2 . ;i )- x 


(T - r2\-1 - 1 (I-Bp-'BlSoV-Bl)- 1 
1 r) 2 1+i tr[(I—B^)~ 1 B^8o] ’ 


(A.22) 


and thus 

1 +tr [(/ + K^K^K^} 


= 1 - 3 tr [( 7 “ B r 


1 B r 5o\ + 4 


ltr \(I-B?) 

- T 

2 


^^ BtSoU-Bp-'BrS 0 ] 

1+itr [(I-B?)-iB?5 0 ] 


Now writing B t 5q as ^ <B> 1 (where V’(A) = Ai(A + r)), A = (I — P 2 ) 1 and P 
checks that 

tr[AP(-i/> <g 1 )A(i/} ® 1)] = tr [AB(ip <g> 1)] tr[A(^ (g 1)]. 


B r 5o, one 
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Using this fact, the above identity yields 

i , 17 T, TSO. ts ^-l z,- 1 l-\tv[(I-Bl)- 1 B r 5o] + \ tr[{I - B?)- 1 B?5 0 \ 

1+tr l(I + K ltl K 2j2 ) A 1i2 J =- 1+ i tr[i i_B?)-iBr6 0 ] - 

Using this and (A.21) in (A.19) we deduce that 

2 (1 - \ tr[(/ + Br)- 1 B r 5 0 ]) 2 


l-itr[(/+B r )- 1 J B r 5o] 
1+1 tr 


Pf [J-K] =det[l -B z r 


1 + 4 tr[(7 — B^)~ 1 B^5q] 


Now let <p e (x) = (^) 1 ^ 2 e l2 /( 2e ) for e > 0, so that 


tr[(I + B r ) i B r 5o] = limtr[(7 + B r ) l B r (ip e <g) 1)]. 

e —>-0 


Since tr[^ £ (8) 1] = 1 we have 


1 — 5 tr[(I + B r ) 1 B r 6 0 \ = i ^ lim tr[(7 + B r ) 1 (tp e (g) 1 )] = ± ± tr[(7 + B r ) 1 <5 0 ]. 


- 1 / 


\—11 


£—^0 


In a similar way we get 


1 + itr [(I-B 2 ) 1 B 2 S 0 ] = 5 + ^tr[(I-S 2 ) 1 <5 0 ] = £ + | tr[(7-5 r ) 1 <5 0 ] + | tr[(7+R r ) 1 <J 0 ]. 


?2\-l 


\-ii 


\-ii 


We deduce from the above identities that 

+i 


Pf[J - K] = det[1 - B 
Now we use the identity 


(1 + tr[(J + B r ) ho]) 


2 + tr[(J - -Br) _1 ^ 0 ] + tr[(J + B r )- 1 6 0 \ 
det[/ — B r \ 


= tr [(I + B r ) 1 h 0 ], 


det[J + B r 

proved in [FS05], and the fact, noted in [LDC12], that the same identity holds (with 
identical proof) if B r is replaced by —B r . The consequence is that 


Pf [J - K] 2 = det[l - Bl 


(l 1 det 

[I-Br}\ 2 

l 1 1 det 

I+B t ] > 

r, , detf I—Br 

, detp+Ffo 

“ 1 det [I+B r 

1 det [I-B r ] 


= det [I — B r 


(A.23) 


Since, by [FS05], det [I — B r \ L 2 ([ 0 ,oo)) = -^GOe(t), this proves (A. 14) up to sign. 

In order to determine the sign the basic idea is to argue by continuity and compare 
the two sides in the limit r -+ oo (in the remainder of the proof we will reintroduce 
the subscript in K r ). In order use continuity we will take advantage of the fact that, 
although K r is not trace class, it is easy to turn its Fredholm Pfaffian into that of a 
Hilbert-Schmidt operator (for the definition see [QR14]). In fact, defining a multiplication 
operator M(^)(x) = 1 f+i)) x ) = (1 + x 2 ), it is not hard to check, using the 

fact that |Ai(x)| < ce~ 2x3/2 / 3 for x > 0, that M J K r M is Hilbert-Schmidt (see Example 
2 in Section 2 of [QR14] for the proof of a similar fact). Thus we may use the notion of 
regularized Pfaffians introduced in Section B.3 (see (B.18)). The idea is to study 

Pf 2 [J - M J K r M] = Pf [e~^ JMlKrM (J + JM T K r MJ)e~^ MTKrMJ ]} 

= e ^- ] ( Kr )i : 2<P] Pf [j + JM T K r MJ}. 


The first equality is by definition of the regularized determinant. For the second one, where 
< i i and 0 _1 denote the corresponding multiplication operators, we are using the fact that 
the operator appearing in the exponent is trace class (see Example 3 in Section 2 of [QR14] 
for the proof of a similar fact) to show that det[e~ 2 MTKrMJ ] equals the exponential; the 
identity follows now from (B.16) (the fact that Pf[J + JM T K r MJ] defines a convergent 
series follows in the same way as for Pf[J — K r ]). 

A similar argument to the one that shows that M T K r M is Hilbert-Schmidt shows that 
M t K r M is continuous in r in Hilbert-Schmidt norm, which implies by (B.19) that the left 
hand side of (A.24) is continuous in r. Moreover, similar arguments show that 4>~ 1 (K r )i : 2<t> 
is continuous in r in trace class norm, and thus the exponential in (A.24) is continuous 
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in r. We deduce that Pf[J + JM T K r MJ] = Pf[J — K r \ is continuous in r and thus, 
since -Fgoe(t) > 0 for all r, this and (A.23) imply that there is a a E {—1,1} such that 
Pf [J — K r \ = <7 -Fgoe (t) for all r. Using again (B.19) and the fact that lim r _>. 0 0 M T K r M = 


M t KqqM in Hilbert-Schmidt norm with /Foo(Ai, A 2 ) = 


0 0 

0 4 sgn(A2—Ai) 


, we get lim^oo Pf [J— 


M J K r M] = Pf [J—M J K^AI] = 1. Since lirn^oo -Fgoe(t) is also 1, this implies that <7 = 1 
and finishes the proof of (A.24). □ 


Appendix B. Fredholm Pfaffians 

The purpose of this section is to provide a brief discussion about the theory of Fredholm 
Pfaffians. These were introduced by Rains in [RaiOO], who stated its main properties. Some 
of the issues discussed in Section B.3 are not directly motivated by what is used in the main 
text, but we chose to include them in the hope that they will help clarifying the notion of 
a Fredholm Pfaffian and its relation to the Fredholm determinant. 


B.l. Pfaffians. Let A be a 2n x 2n skew-symmetric matrix. By basic spectral theory the 
eigenvalues of A come in pairs ±iA (for real A), which then implies that det(A) can be 
written as the square of a polynomial in the entries of A. The Pfaffian of A is defined (up 
to sign) to be this polynomial, and it has the following explicit form: 

1 n 

Pf (^) = 2^! S S n ( a ) II A o(2j-l),a{2j)- (B.l) 

3= 1 

When A has odd size the same arguments imply that A has 0 as an eigenvalue, which leads 
to define Pf(A) = 0. The fact that the definition of the Pfaffian through (B.l) satisfies 

Pf(A) 2 = det(A) (B.2) 

is not at all obvious; we refer the reader to Section 3 of [DG09] which contains three 
different proofs. 

Another basic fact about Pfaffians (a proof of which can also be found in [DG09]) is the 
following: Given two matrices A, B of size k x k, with k even and A skew-symmetric, 

Pf(BAB J ) = det(.B) Pf(A). (B.3) 


A particular case of this identity, which we use repeatedly in Section 4 and below, is the 
following. Let A, B, U be given matrices of size k x k, with A and B skew-symmetric, and 
D 1 , D 2 be diagonal matrices of the same size. Then 


k 

n(£> 1 )o, a (£> 2 )o, 0 pf 

a= 1 




pf 


(D 1 0 \f A 

Vo D 2 )\-U 


U\ (D 1 0 \ 

B) V 0 D 2 ) 


(B.4) 


We collect here some further facts about Pfaffians of skew-symmetric matrices which 
were used in the main text. The first one was used to extend certain formulas from an even 
to an odd number of variables in Section 4.2, while the next two involve certain integration 
identities which were useful in the derivation of the Pfaffian formulas in Sections 4 and 5. 


Lemma B.l. Fix k € Z>o and consider two skew-symmetric k x k matrices A and B, two 
(column) vectors U, V of size k, and define the k x k matrix D = UV J . Then 


Pf 


A 

-D t 


D 

B 



Pf [A] Pf [B] 


Pf 

A U 

Pf 

B V 

-U T 0 

-U T 0 


if k is even 
if k is odd. 


(B.5) 
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Proof. Let us first consider the case k even. More generally, we will consider the case where 
A is k\ x k\, B is k 2 x k 2 and D is k\ x k 2 , with k\ and k 2 even. If A or B are singular 
then clearly both sides vanish, so we will assume that A and B are invertible. We have 

2 


Pf 


A D 
-D t B 


= det 


A D 
-D t B 


= det [A] det [B + D J A~ 1 D}. 


Now D J A~ 1 D = V(U T A~ 1 U)V J , and the middle factor is a scalar which, by the skew- 
symmetry of A _1 , has to be zero. This shows that 

2 


Pf 


A D 
-D t B 


= Pf[A] 2 Pf[B] 2 , 


(B.6) 


which implies the desired identity up to sign. To determine the sign, replace D by eD for 
e > 0. Since Pf _^, T £ g is a continuous function of e which does not vanish for e > 0 

(recall that we are assuming that A are B are not singular) we see by (B.6) that its value is 
either always Pf[M] Pf[B] or always — Pf[v4] Pf[B], Now taking e = 0 the Pfaffian becomes 
Pf [ 0 b] = which shows that (B.6) holds without the squares as desired. 

Now for the case A; odd, consider the matrix 


E = 


By the previous case, Pf [Tv] = Pf 


r a 

D 

U 

U - 

-d t b 

—V 

-V 

-u T v T 

0 

1 

L -U T V T 

-1 

0 . 

d' 

T B 

Pf[ 

0 1 
-1 0 

] = 


A D 

-d t b 


so we may work with 


E. Now permuting the second and third rows and columns of E does not change the value 
of the Pfaffian, and leaves us with 


E' = 


A U D U 
-U J 0 P T 1 
-D J —V B —V 
_-U T -1 V T 0 


Noting that 
follows. 


D u 
-V T 1 


is rank one we see that we are in the first case (k even) and the result 

□ 


Lemma B.2. Consider three skew-symmetric kernels A, B, C defined on some measurable 
space (X,n). Then, assuming that all integrals converge, we have 


-! 

n\ Jx 2 > 


M® 1 2n (dx ) Pf [(A + B)(x a ,x b )] 2 a n b=1 Pf [C(x a ,x b )] 2 a n b=1 


E 


H mn (dx) Pf [A(x a ,x b )]f b=1 


ki,k2>0, k\+k2=n 


(2k 1 )\(2k 2 )\ 


lx 21 


x Pf [B(x a ,x b )]l k ^ +1 Pf [C(x a ,x b )] 2 ™ b=1 . 


Proof. Using several times the definition of the Pfaffian, the fact that the Pfaffian is skew- 
symmetric, and the symmetry of the integration measure, the left hand side of the identity 
equals 

z ln - j Jx <r£S 2n a =1 

1 P ^ . 

- / ; p, s2n (dx) Y\(A + B)(x 2a -i,X 2 a)Pt [C{x a ,x b )] 2 ™ b=1 

lm •* x2n a =1 


2 n n 
1 


2^-7 [ H® 2n {dx) ^ nA{x2a-l,X2a)Y[ B ( X ‘2a-l,X2a)P{[ C ( X °’’ X b)] 2 2 b= l 

dX^ n T ,— ri ~i t ^ t 

1 


2 n n! 


E 

m =0 


lx 


/C{l,...,n} aG/ atfzl 

m n 

/J® 2n (dx) Yl A ( x 2 a-l,X 2 a) B(x 2a -1, X 2a ) Pf [C(x a , X 6 )] ^ =1 


a= 1 


a=m -\-1 
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where we have used the fact that, for |/| = m, the permutation which maps the ordered 
m-tuple formed by the elements in the set U ae/ {2 a — l}U{2a} into the m-tuple (1,..., 2m) 
is even, which implies that the Pfaffian of the kernel C does not change after reordering 
the variables. Using again the symmetry of the integration measure and the antisymmetry 
of the Pfaffian we get that the above 


n .. m 

^ 2 n m\(n - m)\ J x m (>bj! ^ H A ( x M 2 a-i), x ai{2a) ) 

m= 0 v y A \ J (T1 g5 2m a =l 


1 


(2 n — 2m)! 


&l£S2m a=l 
n 

e n (%cr 2 (2a —1) ? *^cr2(2a)) ^§n((Ti) Pf ^fr)] a b = \ 

cr 2^S2(n-m) Q>=m +1 


n 1 r 

E (2m)!(2 „ _ 2m)! f x ,pf 


which gives the desired result. 


x Pf [5(s a , x b )] a " 6=2m+1 Pf [C(x a , ® 6 )] a " 6=1 , 

□ 


Our second integration formula can be regarded as a certain Pfaffian version of the 
Andreief identity [And83] (sometimes referred to as the generalized Cauchy-Binet identity): 

Lemma B.3. Let (A, p) be some measurable space and suppose that for every x € X k 
the matrix [A a ^(x a , x b )] a b=1 is skew-symmetric. Let B be another k x k skew-symmetric 

matrix and for x € X k consider another matrix \U aib (x b )\ k ab=l - Finally consider functions 
<p a defined on X. Then 



= Pf 


k 

(dx) JJ <?!> a (® a )Pf 

a= 1 

[ fx 2 T(dx)h(dx')f> a (x)c/> b (x')A atb (x, x’ 
- [ fx h( dx ) M x ) U b,a( X )\ a,b=l 


[A a ,b(.X a , £&)] ab=1 

[ — U b) a(x b )\ b=l 


[U a ,b(x a )\ ab=1 

[ B -A k a , b =l \ 

)]*6=1 [fxh(dx)(j) a {x) u a)b {x)\ 


k 

a, 6=1 


provided that all integrals converge. 


Proof. Use (B.4) on the left hand side of the claimed identity and expand the Pfaffian using 
its definition to see that it equals 

J k h{dxi) ■ ■ ■ n{dx k ) 

(J ^2k 

| -^cr(2a— l),cr(2a) fa a (2a—1) 1 %cr(2a) )^cr(2a— 1) facr(2a—l) )0cr(2a) faa(2a )) 

a(2a—l)<k 
a ' cr(2a)<k 

^ | U(j(2a— l),cr(2a )—k faa(2a— 1) )4 > cr(2a— 1) facr{2a— 1)) 

a(2a—l)<k 
a ‘ cr(2a)>k 

^ | U a (2a)—k,cr(2a— 1) faa(2a) ) ( / ) cr(2a) fa a" {2a )) | ^cr{2a— 1)— k,cr(2a)—k • 

a(2a—l)>k cr(2a—l)>k 

a ‘ tj{2a)<k a ‘ cr(2a)>k 


Since each x a (a = 1,..., k) appears in one and only one of the above factors, the integrals 
can be brought inside the corresponding factor, and now forming the resulting Pfaffian 
gives the identity. □ 
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B.2. Fredholm determinants. It is instructive to review briefly the theory of Fredholm 
determinants on L 2 spaces before discussing Fredholm Pfaffians. For more details see 
Section 2 of [QR14] . 

Let A by an m x m matrix, write [m] = {1,... ,m } and let (^) denote the family of 
subsets of [m] of size l. A standard calculation shows that, for A 6 C, 

m 

det(7 + AA) = E* E det(Asxs) 

£=0 5e( [ 7 ] ) 


with self-explanatory notation. This suggests a way of extending the determinant to the 
infinite dimensional case. In the case of L 2 spaces it leads to the following. 

Let (X, E,/r) be a measure space and K: X x X —> R be a kernel which defines an 
integral operator acting on L 2 (X) through Kf{x ) = [ x fi(dy) K(x,y)f(y). The Fredholm 
determinant of K (on L 2 {X)) is defined to be the (formal) power series 


°°^ yk r j- 

det 1 1 + A K] t2(y ^ — / y (dx\)... y, (dx n ) det K(x a ,x b ) 

k =0 K ' ^ Xk 


- a,b= 1 


(B.7) 


We will omit the subscript L 2 (X) in the determinant when no confusion can arise. This 
identity can be regarded as a numerical identity whenever the right hand side is absolutely 
convergent. This is the case, for instance, whenever \K(x,y)\ < C for all x,y € X and 
some C > 0, thanks to Hadamard’s bound. 

An interesting class of operators for which the Fredholm determinant power series is 
absolutely convergent is the family of trace class operators. We recall that an operator 
K: L 2 (X) —> L 2 (X) is said to be trace class if it has finite trace norm: 


I All -, : = 


E j x ‘ 

n> 1 JX 


y,{dx) ^ n {x)\K\il: n {x) < oo, 


where (V>n)n>o is any orthonormal basis of L 2 {X) and |AT| = \f~K*K is the unique positive 
square root of the operator K*K. Such an operator is necessarily compact, and in this 
case one has 

det[J + zK] = ]J(i + z\ k ) (B.8) 

fc>i 

where the A^’s are the eigenvalues of K. This identity (known as Lidskii’s Theorem) is 
highly non-trivial, and provides one possible way to extend the definition of the Fredholm 
determinant to trace class operators on a general separable Hilbert space. A nice property 
of the Fredholm determinant restricted to trace class operators is that it is continuous: for 
K\ and A^ trace class one has 

|det[I + Ad] — det[7 + AT 2 ]| < ||Ad - iL 2 ||ie l|Kl||l+l|i ^ 2||l+1 . (B.9) 

This inequality follows from a rather general (and simple) argument based on the analiticity 
in A € C of the function det [I + XK] and the inequality |det[7 + AA']| < e^lll^H 1 (which 
essentially follows from (B.8)), see Theorem 3.4 of [Sim05] or Corollary II.4.2 of [GGK00] 
(see also the discussion following (2.9) in [QR14] ). 

Another useful fact about Fredholm determinants is the so-called cyclic property: if 
K \: A 2 (Xi) —> L 2 {X 2 ) and K\ : L 2 (X 2 ) —> L 2 {X i) then 

det [I + K\K 2 ] = det [7 + K 2 K{\ 


whenever the two sides are absolutely convergent (e.g. if both K\K 2 and K 2 K\ are trace 
class). An important consequence of this fact is the ability to conjugate an operator without 
changing its Fredholm determinant: 

det [7 + V~ X KV\ = det [7 + AT] 


(B.10) 
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whenever both sides make sense. This is often used to replace K by a conjugate trace class 
kernel. For more on all this and on extensions to other separable Hilbert spaces see [Sim05; 
GGKOO], 

Finally let us briefly recall the notion of regularized determinants, since a similar notion 
will be introduced below in the Pfaffian case. Suppose that K is a Hilbert-Schmidt operator 
(see [QR14] for the definition) and recall that the product of two Hilbert-Schmidt operators 
is trace class. We define the regularized determinant of K as 

det 2 [J + K\ = det[(J + K)e~ K \. 

Note that (/ + K)e~ K — I = Yln >2 ~ -AT n is a trace class operator, so the def¬ 

inition makes sense for any Hilbert-Schmidt operator K (one can also define regularized 
determinants of higher order, but we will not do it here). It can be shown that the analog 
of (B.8) for the regularized determinant is 

det 2 [I + K] = J](l + A fc )e- Afc , 

k> 1 

which further implies that |det 2 [F + K\\ < e^ K ^. Using this bound one can obtain the 
analog of (B.9) (by the same argument): if K\ and K-> are Hilbert-Schmidt operators then 

|det 2 [F + K\] — det 2 [/ + Kf[\ < - K 2 \\ 2 e^ K ih+\\ K ih+i)\ (B.ll) 

which in particular gives continuity of the regularized determinant with respect to the 
Hilbert-Schmidt norm. Another way in which this notion can be useful is the following. 
Suppose one knew that both det[I + K] and det[e _A ] are given by asolutely convergent 
series. Then det [(I+K)e~ K ] = det [/©AT] &et[e~ K ] (which follows from the general property 
det[(I + K\) (I + K‘))\ = det [I + K{\ det [I + A 2 ])- The left hand side can be controled in 
terms of the Hilbert-Schmidt norm of K, so if one has some additional control on det[e~ A ] 
(observe, in particular, that if K is nice enough, e.g. trace class, then det[e _A "] = e — tr l^l) 
this can be used to control det [I + K). For much more on this see [Sim05; GGKOO]. 


B.3. Fredholm PfafRans on L 2 spaces. We turn now to Fredholm PfafRans. Given 
n € Z>i we define J to be the 2 n x 2 n block-diagonal matrix which has 2x2 blocks on 
the diagonal, all equal to [ _° 1 q ]. If A is a 2n x 2n skew-symmetric matrix then one can 
show that 

m 

Pf(J + XA) = E a ' e Pf(^5x 5 )- 

£=0 Se( [ 7 ] ) 

This is, of course, the Pfaffian analogue of (B.7), and suggests the following infinite di¬ 
mensional extension [RaiOO] (note that the above sum actually involves only even £, since 
Pf(As'xs) = 0 otherwise) 

Consider a skew-symmetric 2 x 2-matrix kernel 


K{ Ai,A 2 ) 


U,l(^lA2) Kl,2(Ai,A2) 
— A’l,2(A2,Al) A'2,2 (Ai,A2) 


(the skew-symmetry condition in this case translates into AT a ,a(Ai, A 2 ) = —A"a,a(A 2 , Ai) for 
a = 1, 2). We regard K as an integral operator acting on / G L 2 (X ) © L 2 (X) as follows: 


(Kf)b(x) = Y] I v(dy)K b!a (x,y)f a (y) for b = 1,2. 

a=l d X 


For any such kernel we define its Fredholm Pfaffian on L 2 (X) as the (formal) power series 

Pf[,J + \K] l2(x) = YjT [ pidxx)---ii(dx k )P{[K(x a ,x b )] k ab= i ( B - 12 ) 


k> 0 
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whenever the right hand side is convergent. As for Fredholm determinants, we will usually 
omit the subscript L 2 (X) in the Pfaffian 6 . 

Observe that for a 2n x 2 n skew-symmetric matrix A we have by (B.3) that 

Pf (J + A) 2 = Pf(J(J + J~ l A)f = det(J) det(J + J~ l A) = det(J - JA), (B.13) 

where we have used the facts that det(J) = 1 and J 2 = —I. The following result extends 
this to a relationship between Fredholm determinants on L 2 (X) © L 2 (X) and Fredholm 
Pfafhans on L 2 (X). 

Proposition B.4. For any skew-symmetric 2 x 2-matrix kernel K and any A E C, as long 
as the integrals in the series (B.7) and (B.12) are convergent, we have 

Pfb + V© W = det[/- \JK]^ (X) ^ {X) , (B.14) 

where the identity is in the sense of formal power series. 

Proof. Let Pf[J + A K'j = ^ \ n b n and det[/ — A JJF] = ^A n a n . If K is finite-rank 
then (B.13) implies that a n = J2j =o ("j bjb n ~ 3 . In the general case the kernel K can be 
approximated by discretization (as explained in Section 2 of [QR14]) to see that the relation 
between the coefficients of the power series still holds, which proves the identity. □ 


Of course, (B.14) is to be regarded as a numerical identity whenever both sides are 
absolutely convergent. As we will see next, this is the case when K is trace class. 

Proposition B.5. If K is a skew-symmetric 2x 2-matrix kernel which defines a trace class 
operator on L 2 (X) © L 2 (X), then both sides of identity (B.14) define absolutely convergent 
series, and in particular the identity is numerical. Moreover, if the Fredholm Pfaffian series 
of K is given by Pf[J + XK] = En>o a ™A n , then E n >o l a n,A n | < 1,1 . 


Proof. The fact that the right hand side of the identity defines an absolutely convergent 
series follows from our discussion about Fredholm determinants and the fact that \JK\ = 
\K\, so that || JK\\i = ||AT||i. Now suppose that A is a 2n x 2n skew-symmetric matrix and 
let z € C. Since J + zA has pure imaginary eigenvalues in complex conjugate pairs and 
det(J + zA) = — det(7 — zJA), we know that the eigenvalues of JA are real and they all 
have even multiplicity. This means that det (I — zJA) = Pf (J + zA) 2 = ]/[ a (l — z\j) 2 where 
the Aj’s are the eigenvalues of JA (appearing with half their multiplicity). This identity 
can be extended by approximation to K since K is trace class, and this determines the 
Fredholm Pfaffian series for K up to sign: 


Pf[J + zK] = ^2 a nZ n = (J JJ(1 ~ z\j) 

n> 0 j 

for a € {—1,1}. In particular, this means that the n-th term of the Fredholm Pfaffian 


series is given by a n = az n ^ 




Aji Now 


E 


h<-jn A A ''' Aj„ 


< U\\\J K \\Ti 


by Lemma 3.3 and (3.8) in [Sim05] (the 1/2 is because each eigenvalue is counted with half 
its multiplicity). Since || JK\\\ = ||77||i as explained above, we have proved that 


El 

n> 0 


a n z 


7l>0 


1 




n (1 | 

2 I 


K 


Ljra _ e \\A\\ K h' 


□ 


An important consequence of the last result is that the Fredholm Pfaffian restricted to 
trace class operators is continuous: if K\ and K 2 are trace class skew-symmetric 2x2- 
matrix kernels then 

|Pf[J + K\] — Pf[J + AT 2 ]| < ||Ad -R 2 ||ie^ (l|A ' l||l+ll ^ 2||l+1) . (B.15) 

^Observe that, although K acts on L 2 (X) © L 2 (X), we are declaring this to be the Fredholm Pfaffian 
on L 2 (X). This is just a matter of convention, so we make this choice for notational convenience. 
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Like (B.ll), this inequality can be proved in the same way as (B.9), based now on the 
inequality |Pf[J + AT]| < implicit in Proposition B.5. 

We mention also two Fredholm Pfaffian analogs of (B.10). The first one follows from 
(B.3) by approximation and says that 

Pf[(I + L J ){I + K)(I + L)\ = det[J + L\ Pf[J + K] (B.16) 

as long as both sides define absolutely convergent power series. For the second one, let 
M be a symplectic 2x2 matrix kernel. Then, as long as both sides define absolutely 
convergent power series, we have 

Pf[J + A K\ = Pf [J + A M J KM\. (B.17) 

The proof of (B.17) is elementary: 

Pf[J + XM J KM] = ^2 JJ [ fi(dx 1 )---fi(dx k )Pf[(M J KM)(x a ,x b )] k ab=1 
k >o K ' J xk 

= E4r/ »(dxi) ■ ■ ■ n(dx k ) det[M(x a ,x b )] k ab=1 Pf[K(x a ,x b )] k ab=1 = Pf [J - K], 
k >o Jxk 

where we have used (B.3) and the fact that the determinant of any symplectic matrix 
equals 1. 

One can in principle use (B.17) to replace K by a trace class kernel and then use the 
(better-developed, and in a sense simpler) theory of Fredholm determinants to study the 
Fredholm Pfaffian of K. On the other hand, even if this is not possible one can still use 
Proposition B.4 to reduce the study of a Fredholm Pfaffian to a Fredholm determinant, as 
long as they are both convergent. This is in fact the situation we are in in Section A.3, 
where we are dealing with a kernel K for which there is no obvious symplectic 2 x 2-matrix 
kernel such that M T KM is trace class, but the reduction to a Fredholm determinant is 
still very useful. 

Let us finish by introducing a notion of regularized Fredholm Pfaffians. Suppose that K 
is a skew-symmetric 2 x 2-matrix kernel which defines a Hilbert-Schmidt operator. Observe 
that J is symplectic, and thus by (B.17) we have 

Pf[J + K] = Pf[J + J t KJ] = Pf [J - JKJ}. 

Observe also that ( JK) T = KJ. In view of these two facts we define the regularized 
Pfaffian of K as 

Pf 2 [J + K] = Pf[(J - JKJ)e KJ ] = Pf[e^ JK (J - JKJ)e^ KJ }. (B.18) 

The second equality actually holds at the level of the operators inside the two Fredholm 
Pfaffians (as can be checked easily using the series expansion of the exponential) and has 
the advantage of making the skew-symmetry of the argument more apparent. A calculation 
shows that (J—JKJ)e KJ — J = — Yln >2 and thus the right hand side of (B.18) 

makes sense as the Fredholm Pfaffian of a trace class operator (since JK is Hilbert-Schmidt, 
and thus ( JK) n is trace class for any n > 2). Note also that 

Pf 2 [J + K] 2 = det 2 [/- JK}. 

Additionally, if K\ and A" 2 are Hilbert-Schmidt then 

|Pf 2 [J + K{\ — Pf 2 [J + Kf\ \ < ||A'i - A 2 || 2 e3 ( l |A ' l||2+l|A ' l||2+1)2 (B.19) 

(with the same proof as (B.9) , now based on the inequality | Pf 2 [ J+K] \ = y^det 2 [/ — JK] \ < 



Observe that if one knew in addition that K is trace class, then by (B.16) one would 
have 

Pf 2 [J + K] = det[e^ JA ] Pf[J - JKJ] = e ^ tr[JK] Pf[J + K). (B.20) 
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As in the case of the regularized determinants discussed in Section B.2, one situation in 
which the introduction of regularized Pfaffians can be useful is when K is Hilbert-Schmidt 
but not trace class, but on the other hand some additional control on K is available which 
makes it possible to show that the first or both equalities in (B.20) hold. In such a case, 
since the left hand side can be controlled by the Hilbert-Schmidt norm of K (which is 
the same as that of JK ), this identity provides a possible route for controlling Pf[J + K]. 
This idea was used at the end of the proof of Proposition A.4 to use continuity in r of 
a certain kernel K r with respect to the Hilbert-Schmidt norm, in order to upgrade the 
identity Pf [J — I\ r ] 2 = det [I + JK r \ to Pf [J — K r \ = \J det[I + JK r \. 
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